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ABSTRACT 

Context. Nitrogen is a key element to test the impact of rotational mixing on evolutionary models of massive stars. Recent studies 
of the nitrogen surface abundance in B-type stars within the VLT-FLAMES survey of massive stars have challenged part of the 
corresponding predictions. To obtain a more complete picture of massive star evolution, and to allow for further constraints, these 
studies need to be extended to 0-stars. 

Aims. In this paper, we investigate the N iv/i4058 emission line formation, provide nitrogen abundances for a substantial 0-star sample 
in the Large Magellanic Cloud, and compare our (preliminary) findings with recent predictions from stellar evolutionary models. 
Mefhods. Stellar and wind parameters of our sample stars are determined by line profile fitting of hydrogen, helium and nitrogen 
lines, exploiting the corresponding ionization equilibria. Synthetic spectra are calculated by means of the NLTE atmosphere/spectrum 
synthesis code fastwind, using a new nitrogen model atom. We derive nitrogen abundances for 20 O- and 5 B-stars, by analyzing all 
nitrogen lines (from different ionization stages) present in the available optical spectra. 

Results. The dominating process responsible for emission at Niv/14058 in 0-stars is the strong depopulation of the lower level of the 
transition, which increases as a function of M. Unlike the N iii triplet emission, resonance lines do not play a role for typical mass-loss 
rates and below. We find (almost) no problem in fitting the nitrogen lines, in particular the 'f ' features. Only for some objects, where 
lines from N iii/N iv/N v are visible in parallel, we need to opt for a compromise solution. 

For five objects in the early B-/late O-star domain which have been previously analyzed by different methods and model atmospheres, 
we derive consistent nitrogen abundances. The bulk of our sample 0-stars seems to be strongly nitrogen-enriched, and a clear correla- 
tion of nitrogen and helium enrichment is found. By comparing the nitrogen abundances as a function of v sin ; ('Hunter-plot') with 
tailored evolutionary calculations, we identify a considerable number of highly enriched objects at low rotation. 
Conclusions. Our findings seem to support the basic outcome of previous B-star studies within the VLT-FLAMES survey. Due to the 
low initial abundance, the detection of strong Nitrogen enrichment in the bulk of O-stars indicates that efficient mixing takes place 
already during the very early phases of stellar evolution of LMC O-stars. For tighter constraints, however, upcoming results from the 
VLT-FLAMES Tarantula survey need to be waited for, comprising a much larger number of O-stars that will be analyzed based on 
similar methods as presented here. 

Key words, stars: early-type - stars: abundances - line: formation - stars: atmospheres - stars: winds, outflows 

1 . Introduction e.g., lHerreroll200l iHerrero & Lennon ''2004' and 'Moref '2009|). 

Several recent studies based on the nitrogen diagnostics per- 

One of the main products of the VLT-FLAMES survey of mas- ^^^^^^ ^^^^^^ FLAMES I survey have severely challenged 

sive star^ (hereafter FLAMES I) was the determmation of hght ^^e predicted effects though. In this context, nitrogen is a key el- 

element abundances from statistically sigmficant samples of ement to test the predictions of rotational mixing, since it should 

Galactic, Lai-ge and Small Magellanic Cloud (LMC, SMC) B- ^^^^^^ ^^^.^^^^^ enriched at the stellar surface of rapidly ro- 

^ stars, covering a broad range of rotational velocities. j^jj„g ^j^^, ^^^^^^ j„ ^^^^^ g^^jy evolutionary phases, whilst 

The inclusion of rotational mixing mto massive star evo- for slow rotation almost no enhancement should occur before 

lution (e.g., iHeger & Langeij |2000t |Meynet & Maedeij |2000t the red supergiant phase. Actually, a significant number of 

iBrottet al.||2011a[) brought better agreement with spectroscopic both un-enriched fast rotators and highly enriched slow rotators 

analyses that provide evidence for a characteristic enrichment have been foun d within the population of LMC core-hydrog en 

of helium and nitrogen in many eai'Iy-type stai's (reviewed by, burning objects (iHunter et alJl20()la 120091: iBrott et aI.ll2()TTM) FI 

These results imply that standard rotational mixing might be not 

* Based on observations collected at the European Southern dominant, and/or that other enrichment processes might be deci- 

Observatory Very Large Telescope, under programmes 68.D-0369, sive as well (Brott et al 201 iB) 

171.D-0237 (FLAMES) and 67.D-0238, 70.D-0164, 074.D-0109 rr, . . • ' u j- 

(UVES) further constrain these nndings and to provide a general 

T3 J 1 1 ui ■ 1 f ■ f f picture of massive star evolution, these studies need to be ex- 

Appendix A, B and C are only available m electronic torm at ^ ' 

|http://www.edpsciences.org 
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see Evans et al] 120061 for an introductory publication and ^ as well as slowl y rotating, highly enriched supergiants, discussed 

lEvans et al..i2008i for a brief summary on the outcome of this project bv lVink et alj ( 120101) 
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tended to O-type stars. Because of their shorter Hfetimes, the 
time-range where this enrichment takes place can be narrowed 
down, and one might be able to constrain the mixing scenario 
even better than it is possible from B-stars alone. In this re- 
spect, the LMC is an ideal testbed, since the nitrogen baseline 
abundance is low and even a strong enrichment is easier to mea- 
sure/confirm than, e.g., in the Milky way. 

Interestingly, most previous abundance studies of massive 
stars are strongly biased towards intermediate and early type 
B-stars. Indeed, when inspecting the available literature, metal- 
lic abundances, in particular of nitrogen, are scarcely found for 
O-stars. The situation for LMC objects is even worse, and data 
for only a few supergiant s (Pauldrach e t al.ll994l:I Crowther et aP 
l2002t lEvans et al.ll2004 and giants (Walborn et al.l |2004 are 
available. One of the reasons for this lack of information is that 
the determination of nitrogen abundances is a non-trivial task, 
due to the complexity of N m/N iv line formation related to the 
impact of various processes that are absent or negligible at cooler 
spectral types where N ii is the dominant ion. 

To provide more insight into this matter, we started a se- 
ries of publications dealing with nitrog en spectroscopy in O- 
type s tars. In the first paper of this series (iRivero Gonzalez et al.l 
I2OI iL hereafter Paper I), we concentrated on the forma- 
tion of the optical Niii emission lines at AA4634 - 4640 - 
4642, which are fundamental for the definition of the dif- 
ferent morphological 'f -classes (Walbor n 1971b). It turned 
out that the c anonical explanation in te rms of dielectronic 
recombination dMihalas & Hummed [T973h no longer or only 
partly applies when modern atmosphere codes including line- 
blocking/blanketing and winds are used. The key role is now 
played by the stellar wind, which induces a (relative) overpopu- 
lation of the upper level of the transition, via pumping from the 
ground state rather than by dielectronic recombination as long 
as the wind-strength is large enough to enable a significantly ac- 
celerating velocity field already in the photospheric formation 
region. 

The main goal of the present paper is to provide nitro- 
gen abundances for a considerable number of O-stars in the 
LMC. For this purpos e, we use the corresponding sample from 
iMokiem et al.l (I2007al hereafter Mok07), mostly based on ob- 
servations within the FLAMES I survey. So far, this is the 
largest sample of O-stars studied in the LMC by means of 
quantitative hydrogen and helium line spectroscopy, and allows 
us to determine nitrogen abundances for a significant number 
of objects. However, its size is still not comparable with the 
amount of corresponding B-stars, and does not allow us to ex- 
tend the B-star results (that challenged rotational mixing) to- 
wards the O-star domain in a statistically sufficient way. Rather, 
it will yield a first impression on potential problems. A statisti- 
cally significant analysis wil l become possible within the VLT- 
FLAMES Tarantula survey (lEvans et al.ll2011l 'FLAMES IF), 
which provides an unprecedented sample of 'normal' O-stars 
and emission-line stars. 

This paper is organized as follows. In Sect. |2] we describe 
the tools used to determine nitrogen abundances, both the at- 
mospheric code and the nitrogen model atom. In particular, we 
study the formation of the Niv/14058 emission line in parallel 
with the Niv/i6380 absorption line. Sect. [3]presents the stellar 
sample and the observations used within this study. The proce- 
dure to determine stellar and wind parameters together with ni- 
trogen abundances is outlined in Sect.|4] In Sect.|5] we comment 
in detail on the individual objects. Sect. |6]provides a discussion 
of our results, and Sect. Q summarizes our findings and conclu- 



2. Prerequisites for nitrogen diagnostics 

2.1. The code 

For this work, we use a recently updated version (vlO.l) 

of the atmo sphere/line formation code fastwind (see 

ISantolava-Rev et al. 1997 and Puis et alj |2005| for previous 
versions). This code has been specifically designed for the 
optical and IR spectroscopic analysis of hot stars of spectral 
types early A to O, and accounts for NLTE conditions, spherical 
symmetry and mass-loss. The current version incorporates a 
variety of updates and improvements compared with previous 
versions, which are briefly summarized in the following. 

At first note that fastwind differentiates between so-called 
'explicit' and 'background' elements, where the former are those 
used as diagnostic tools (in the present context: H, He and N) and 
are treated with high precision, by detailed atomic models and 
by means of comoving frame transport for the line transitions. 
The background elements (i.e., the rest) are used 'only' for the 
line-blocking/blanketing calculations, and have been treated so 
far by means of the Sobolev approximation. Though this is rea- 
sonable in the wind regime, the Sobolev approximation becomes 
doubtful in regions where the velocity field is strongly curved, 
which is the case in the transition zone between photosphere 
and wind. As we have convinced ourselves, the induced errors 
are not important for the background radiation field, but they can 
have a certain influence on the temperature structure. Applying 
the Sobolev approximation in regions with a strong velocity field 
curvature results, on average, in too strongly populated up per 
levels of line transitions (see, e.g.. Santolava-Rev et al.l ll997h . In 
turn, this leads to overestimated heating rates, which can result 
in too high temperatures in the transition region (and sometimes 
even below). To avoid this problem, the new fastwind version 
treats also the most important lines from the background ele- 
ments in the comoving frame. 

A second modification refers to the photospheric line ac- 
celeration. So far, this quantity (being important for the pho- 
tospheric density stratification - higher line acceleration, lower 
density) has been calculated from the Rosseland opacities, which 
is strictly justified only at large optical depths. In the new ver- 
sion, an additional iteration cycle to calculate the photospheric 
structure is performed, now by using the flux-weighted mean 
from the current NLTE opacities. 

By calculating a large grid of OB-star models, and compar- 
ing with solutions from the previous fastwind version, it turned 
out that both improvements affect mostly dwarfs/giants in the 
effective temperature range 30 kK < Teft < 35 kK. In par- 
ticular the optical Hen lines become stronger, mostly due to 
the somewhat lower (electron-) densities. Interestingly, this is 
just the domain where previous fastw ind solutions showed the 
largest deviations from other codes (ISimon-Dfaz & Stasiriskal 
2008). Comparing the new structures to results from tlusty, ex- 
cellent agreement has been found. For dwarfs/giants with effec- 
tive temperatures outside the 'problematic' region, and for all 
supergiants, the differences to previous results from earlier fast- 
wind versions are small. 

The last major improvement concerns the implementation 
of dielectronic recombination, both for the background and the 
explicit elements, and was already described and successfully 
tested in Paper I. 
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2.2. The nitrogen model atom 

To perform our analysis, we implemented a new nitrogen model 
atom into fastwind, consisting of the ionization stages Nii to 
N V, which has already been used for the calculations performed 
in Paper 1. The level structures of both Ni v and Nv have 
been taken from the WM-basic atomic database (iPauldrach et alJ 
11994 . 



N II has been adapted to fastwind from a previous model ion, 
developed by Becker & Butler (1989). We had no intention 
to develop a 'perfect' Nii model, since most of our analy- 
ses deal with O-star spectra where Nn becomes invisible, and 
particularly because N. Przybilla and co-workers have akeady 
constructed such a mo del (based on an earlier version, see 
iPrzvbilla & ButleJ 12001 ) which will be incorporated into our 
code after release. However, a series of tests were performed to 
ensure the goodness of the model ion, see Sect l2.2.11 

For Nil, we consider 50 LS-coupled terms, up to principal 
quantum number n = 4 and angular momentum / = 3, where all 
fine-structure sub-levels have been packed. Detailed information 
about the selected levels is provided in Table [ATI and Fig. lA.ll 
We account for some hundred permitted electric dipole radiative 
transitions. For the bulk of the transitions, oscillator strengths ar e 
taken from calculations performed bv iBecker & Butleil ( 1989h . 
but for some transitions related to strong Nii lines oscillator 
strengths have been taken from NIST0 Radiative intercombi- 
nations are neglected. Roughly one thousand collisional bound- 
bound tr ansitions are con sidered, with corresponding rates us- 
ing the Ivan Regemorterl ([l962) approximation in the radia- 
tively p ermitted case and following the semi-empirical expres- 
sion bv lAllenl(ll973h in the forbidden one. Radiative ionization 
cross sections have been derived by Becker & Butler' ('1989'), 
and adapted to the representation suggested by Seaton (1958). 
Collisional io nization cross-sections are calculated using the 
ISeatonI ([19621) formula in terms of the photoionization cross- 
section at threshold. 



N III has been already described in Paper I. In brief, it consists 
of 41 packed terms up to n = 6 and / = 4 (doublet and quartet 
system). 

Niv. Also this model ion consists of 50 LS-coupled terms, 
up to principal quantum number n - 6 and angular momen- 
tum I - 4, with all fine-structure sub-levels packed into one 
term. Table IA.2I provides detailed information about the se- 
lected levels. Two spin systems (singlet and triplet) are treated 
simultaneously (Fig. IA.2l i. All allowed electric dipole radia- 
tive transitions between the 50 levels are considered, as well 
as radiative intercombinations, with a total of 520 transitions. 
Corresponding oscillator strengths have been drawn from either 
NIST when available or otherwise from the WM-basic database^ 
Furthermore, we consider roughly one thousand bound-bound 
collisional transitions between all levels, with effective collision 
strengths among the 12 lowest LS-states from R-matrix com- 



firstly described in 



^ http://www.nist.g ov/physlab/data/asd.cfm, 
iKelle her etafllgM) 

see Pauldrach e t alj (11994). In brief, the atomic structure code su- 
perstructure (Eiss ner & Nuss baumer 1969: Eissner 1991) has been 
used to calculate all bound state energies in LS and intermediate cou- 
pling as well as related atomic data, particularly oscillator strengths in- 
cluding those for stabilizing transitions. 



putations bv iRamsbottom et al.l (1 19941) . Transitions without de- 
tailed data and collisional ionizations are treated as in N ii. 

Photoion ization cros s -secti ons have been taken from cal- 
culations by iTullv et al.l (Il990h. via TOPbaseB th e OPACITY 
Project on-line database dCunto & Mendozalll992h . For excited 
levels with no OPACITY Project data available (5g 'G, 5g ^G, 
6s 'S, and 6g ^G, see Table lA!2l l. resonance-free cross-sections 
are used, provided in terms of the Seaton ( 1958) approximation 
with parameters from the WM-basic database. Finally, the most 
important dielectronic recombination and reverse ionization pro- 
cesses are implicitly accounted for by means of exploiting the 
OPACITY Project photo cross-sections. Only for the few lev- 
els with no such data available, we apply the 'explicit' method, 
using the stabilizing transitions on top of resonance-free photo 
cross-section (see, e.g.. Paper I), with corresponding data from 
WM-basic. 



N V. Our model of this lithium-like ion (one doublet spin sys- 
tem) consists of 27 levels, including LS-coupled and packed 
terms up to n = 7 and / = 6 (see Table IA.3I for details and 
Fig. lA.3l for a Grotrian diagram). All allowed electric dipole ra- 
diative transitions are accounted for, with a total number of 102 
radiative bound-bound transitions and oscillator strengths from 
WM-basic. Collisional excitations and ionizations are treated 
as in N ii, whilst photo cross-sections (in terms of the ISeatonI 
'1958 approximation) have been taken from the WM-basic atomic 
database. 

Our complete nitrogen model atom (from N ii to N v) comprises 
178 LS-coupled levels, with more than 1 100 radiative and more 
than 3800 collisional bound-bound transitions. To calculate the 
final synthetic profiles, Voigt profiles are adopted, with central 
wavelengths according to NIST, radiative damping parameters 
from the Kurucz database^ and collisional damping par ameters 

broad ening by electron impact) computed according to lCowlevI 

1971h . 

2.2.1 . Testing tine N ii model ion 

To test our somewhat simple Nii model ionQ we com- 
pared s ynthetic line profiles with corresponding ones from 
TLUSTY (Hubenvl 1998) and detail/surface (Giddings 119811: 
'But ler & Gid'dings[fl985 ). the latter based on the newly devel- 
oped N n model by Przybilla et al. (see Sect. I2.2l i. 

A summary of the various tests is provided in Table [T] We 
started by c ompari ng with line profiles from the BSTAR2006 
model grid (iLanz & Hubenv 2007). This grid has been calcu- 
lated using the model atmosphere code tlusty (IHubenvlll988i : 
iHubenv & Land 19951) . a code that assumes plane-parallel geom- 
etry, hydrostatic and radiative equilibrium, and calculates line- 
blanketed NLTE model atmospheres and corresponding syn- 
thetic profiles. Due to its restrictions, only objects with negli- 
gible winds can be analyzed. 

We used this grid of models as reference because it covers 
the parameter range at which N ii/N iii are the dominant ioniza- 
tion stages, and the N ii lines are clearly visible in the synthetic 
spectra. In addition, the BSTAR2006 grid has been used in nu- 
merous studies aiming at the de termination of ste llar abundances 
and parameters in B-stars (e.g.. lLanz et al.ll2b08l) . 



' http://cdsweb.u-strasbg.fr/topbase/topbase.html 
^ www.pmp.uni-hannover.de/cgi-bin/ssi/test/kurucz/sekur.html 
' important, e.g., for our comparison with B-star nitrogen abundances 
from alternative analyses, see Sect. 16.21 
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Table 1. Nii test series (see text). 



Series 


Atomic model 


Photospheric 


NLTE/line formation 






stratification 




FW 


see Sect. 12.21 


FASTWIND 


FASTWIND 


FW2 


see Sect.l2.2l 


TLUSTY 


FASTWIND 


TL 


TLUSTY 


TLUSTY 


TLUSTY 


Prz 


Przybilla et al. 


KURUCZ 


detail/surface 



A grid of FASTWIND models has been calculated, covering the 
temperature range 20 kK < Teff < 30 kK, using a typical step 
size of 2.5 kK, and gravities representative for dwarfs and gi- 
ants. Since tlusty does not account for the presence of a wind, 
we used negligible mass-loss rates, M- 10"^. . . 10""^ Moyr"'. 
Because fastwind allows for employing external photospheric 
structures, we created a second grid using the tlusty pho- 
tospheric structure, smoothly connected to the wind structure 
as calculated by fastwind, with the usual p velocity law. For 
consistency with the tlusty grid, all models have been cal- 
culated with the 'older' s olar nitrogen abundance, [N] = 7.92 
jGrevesse & Sauva l''1998'). where [N] = log N/H + 12 with re- 
spect to particle numbers. Note that all tests have been performed 
with the complete nitrogen model atom involving the ions N ii to 
N V as described in Sect l2.2l 

It turned out that there are large differences between the syn- 
thetic N II line profiles calculated by tlusty (TL) and fastwind 
(FW), even though the hydrogen/helium lines agree very well 
(except for the forbidden component of He i which is stronger 
in FW-models). The latter code predicts stronger Nii profiles, 
which is also true for the fastwind results based on the tlusty 
photospheric structure (FW2, see Figs. lB.2l to lB.7b . 

For dwarfs, there are almost no differences between the pro- 
files from the FW and FW2 models. This is readily understood, 
since the photospheric stratification of electron temperature, T^, 
and electron density, «e, are essentially the same (see Fig. IB. II 
panel 1, 3, 5), i.e., fastwind and tlusty predict the same struc- 
tures. 

On the other hand, models for giants at higher Teff display 
(mostly) weak differences. E.g., at Teft - 27.5 kK and log ^ = 3.0, 
there is a small disagreement of the electron densities in photo- 
spheric regions (tross < 10"^), even though the temperatures 
agree quite well (Fig. lB.ll panel 6). In particular, the tlusty and 
thus the FW2 structure shows a lower electron density at opti- 
cal depths where the photospheric lines are formed, because of 
a higher photospheric radiative line pressure in this model. 

Because of the lower electron density, the lower recombina- 
tion rates (at T^s - 27.5 kK, Niii is the dominant ion) lead to 
somewhat weaker N ii profiles in the FW2 models compared to 
the FW ones, see Fig. IB. 71 Nevertheless, differences to the pro- 
files as predicted by TL itself are still large, and we conclude that 
the photospheric structure is not the origin of the discrepancies. 

As an independent check, we compared our results with 
spectra calculated by N. Przybilla (priv. comm.) for two of our 
grid models, denoted by Prz in the following. These spectra are 
based on the NLTE/line formation code detail/surface, and the 
Nil model ion recently developed by Przybilla et al.. We con- 
sider this atomic model as the superior one in the present context, 
because large effort on improving and testing the atomic data has 
been spent, and corresponding synthetic spectra perfectly match 
high res olution/high signal-t o-noise observations from various 
B-stars dPrzvbiUa et al.i2008l) . 



Figures IB. 21 and IB .41 show that the agreement between the 
FW and Prz profiles is excellent, which leaves us with the con- 
clusion that there might be problems in the Nii atomic data 
used in the BSTAR2006 grid. Important studies using tlusty 
have been carried out duri n g the past years (e. g., iDufton et al.l 
120061; iTrundle et all 120071; iHunter et al.l l2009l) . including the 
determination of Nii abundances in LMC and SMC B-stars. 
Actually, thes e studies utiliz e d a di fferent model atom, devel- 
oped by AUe nde Prieto et aTl ( l2003h . which has been tested by 
N. Przybilla at our request, with a positive outcome. Thus, the 
aforementioned analyses should be free from uncertainties re- 
lated to a potentially insufficient atomic model. 

2.3. Diagnostic nitrogen lines in fhe optical 

Table |2] presents a set of 5 1 nitrogen lines visible in the opti- 
cal (and adjacent) spectra of OB-stars, along with the position 
of potential blends. Included are the connected levels (for cor- 
responding term designations, see Appendix A) and multiplet 
numbers, to provide an impression of how much independent 
lines are present. 

Lines from N ii (visible in the spectra of B and late O-stars) 
have been selected after careful comparisons with profiles cal- 
culated by N. Przybilla (priv. comm., see Sect. l2.2.TI ). Only one 
of the suggested lines, Nii /13995, is completely isolated, and 
remains uncontaminated even at high rotation rates. Moreover, 
this is one of the strongest N ii lines located in the optical re- 
gion, making it a good choice for deriving nitrogen abundances. 
Other useful lines are N ii/l5667 and /15679, where the former is 
moderately strong and the latter has roughly the same strength 
asNii/}3995. 

The subset comprising lines from Niii has been discussed 
in Paper I. Prominent lines from N iii, N iv and N v are among 
the most well-known features in O-stars, and can be used to 
infer nitrogen abundances as well as effective temperatures for 
the earliest subtypes, from the reaction of the N iv/N v ioniza- 
tionequilibriumQ In a similar line of reasoning, Walborn et al] 
(I2002h used the N iv/l4058 emission line in combination with 
N m/1/14634-464 0-4642 to split the degenerate 03 spectral type 
(IWalbornll971al) into three different types 02, 03, and 03.5, re- 
lying on the Niv/Niii emission line ratio0 Thus, a detailed un- 
derstanding and modelling of N iv/l4058 (together with the N m 
triplet, see Paper I) is mandatory, to safely employ this powerful 
diagnostics. 

NiV/14058 and NiV/i6380 connect the 'neighbouring' levels 
of the singlet series Is^ 2s31 with 1 = s,p,d (levels #8, 9, and 12 
in Table lA. 2b . N IV/14058 (if present) is observed in emission in 
the majority of stars, and has been suggested to be formed by 
photospheric NLTE processes (see below) rather than by emis- 
sion in an extended atmosphere, in analogy to the Niii triplet 
emission. Note that there is no detailed a nalysis of t he line for- 
mation process. So far, onlv lTaresch et al.l (QWTi) and lHeap et al.l 
(I2OO6I) simulated the behavior of this line as a function of ef- 
fective temperature. Heap et al. found emission for this line at 
Teff > 40, 000 K for log g = 4.0 and [N] = 7.92, using the plane- 



^ At the earliest 0-types the standard technique for deriving effective 
temperatures based on Hei and Hen line-strengths becomes difficult or 
even impossible, due to vanishing Hei and rather insensitive Hen lines 
from Teff = 45 kK on. 

^ This c lassification scheme has been criticized by Mass ey et al.l 
( 120041 1200 5^ who found that for stars with similar effective tempera- 
ture and surface gravity the N iv/N m emission line ratio can vary over 
the full range as defined for 02 and 03.5. 
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Table 2. Diagnostic nitrogen lines in the optical (and adjacent) spectrum of early B- and O-type stars, together with potential blends. 



Ion 


Transition 


Multiplet 


Line 


Wavelength(A) 


Blends 


Used 


Nil 


N29 - 


■ N216 


12 


1 


3994.99 


_ 


+ 




N211 


- N219 


15 


2 


4447.03 


On ^4446.81, 4447.67, 4448.19, Oiii ^4447.79 


+ 




N28 - 


■ N215 


5 


3 


4601.47 


On A4602.13, 4603.23, Niii A4604.18 


+ 




N28 - 


■ N215 


5 


4 


4607.16 


On ^4609.44, 4610.20, Niii ^4605.16, Neii ^4606.70 


+ 




N28 ■ 


■ N215 


5 


5 


4621.39 


On ^4621.27, Niii ^4621.04, 4623.05, Sin ^4621.72 


+ 




N28- 


■N215 


5 


6 


4630.54 


Oiii ^4630.77, Niii ^4630.61, Siiv /14631.24 


+ 




N28 ■ 


■N215 


5 


7 


4643.09 


O II A464 1.81, 4643 . 89, N iii A464 1.85 


- 




N212 


- N218 


19 


8 


5005.15 


Oiii/i5006.84 


- 




N214 


- N221 


24 


9 


5007.33 


On: .15006.84 


- 




N28 ■ 


■N214 


4 


10 


5045.10 


Nii/15046.53 


- 




N28 ■ 


■N212 


3 


11 


5666.63 
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Notes. Line numbers for Nn, Niv and Nv refer to important transitions as indicated in Figs. lA.l1lA.2l and lA.3l respectively, and for Nm to the 
corresponding figure (B.l) in Paper I. Multiplet numbers for Nn/Nm are from Moore ( 1975), and for Niv/N v from Moore ( 1971). Lines used 
within the present work are labeled by ' + '. Lower and upper levels of the transitions are denoted by a combination of ion and level number 
according to Tables lA.lilA.2l and lA3] e.g., 'N29' means level #9 of Nn 

parallel atmospheric code tlusty, supporting the idea that the 
emission is of photospheric origin and that velocity fields are 
not required to explain the basic effect. Other arguments for the 
photospheric origin of Niv/14058 are the agreement with other, 
absorption line profiles as a function of v sin unshifted ra- 
dial velocites, and lack of P Cygni profiles (N. Walborn, priv. 
comm.). 



Interestingly, Niv.16380 appears clearly in absorption in O- 
star spectra, and seems to play a similar role as the N iii absorp- 
tion lines at .1.14097 - 4103 in the Nin emission line problem 
(see Paper I). We note that N iv.16380 can be significantly af- 
fected by the presence of two Diffuse Interstellar Bands (D IBs) 
at .1.16376.08, 6379.32 (lHerbigll975l:lKrelowski et an[T995h . the 
latter being stronger if reddening is important. Fortunately, the 
stars analysed in this paper are subject to low reddening, and 
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these DIBs only minorly affect some of the observed N iV/16380 
hnes (e.g. Nll-038,Fig.|C8]l. 

As for the previous pair of Niv lines, NiV/16380 and Nrv 
/14058, also the N iv multiplets around 3480 A and 7103-7129 A 
are formed between levels of the same series (here within the 
triplet system - levels #7, 10, and 1 1 in Table IA.2] |. and seem to 
mimic the behavior of these lines: at least in the earliest O-star 
regime they are prominent features, where the former multiplet 
appears in absorption and the latter in emission. Both line com- 
plexes are widely used in WR-star analyses, and the emission 
in the latter multiplet is a strong feature in most WR spectra. 
The lack of emission at this multiplet and also at N iv/i4058 has 
been used for classification p urposes in different WR-star studies 
(e.g.. lNegueruela & Clarkl20 05). Likewise, the multiplet around 
3480 A has been used bv'Walbo rn et al1(l2004 to infer both T^^ 
and nitrogen abundances for a set of 02 stars. 

The remaining N iv Unes listed, N ivii5200 - 5204 - 5205, 
belong also to the triplet system, and appear, if present, in ab- 
sorption (for O-stars). Unfortunately, most spectra used within 
this study do not cover this spectral region. For the few field stars 
where this range is available to us (see Sect. 13.2b . this multiplet 
is not visible. 

Finally, N v lines at /1/14603-4619 are produced by transitions be- 
tween the fine-structure components of 3s and 3p ^P" (levels 
#3 and 4 in Table lA.3"l l. These doublet lines are strong absorption 
features in the earliest O-stars, showing sometimes extended ab- 
sorption in their blue wings or even pronounced P-Cygni pro- 
files (e.g., Nl 1-031, Fig. lC.7b . revealing that they can be formed 
in the wind. Thus, it is clear that mass-loss and wind-clumping 
will influence the formation of these lines. Besides, we also list 
N V/1/14943 - 4945 which become important, (almost) isolated 
diagnostic lines in the spectra of very early, nitrogen rich O- 
and WNL-stars. Unfortunately, the corresponding wavelength 
range has not been observed for the bulk of our sample stars 
(see Sect 13.2b . whilst no features are visible in the few early- 
type spectra (field stars) where this range is available. Note that 
in order to use these lines, we would need to extend our Nv 
model ion, including high-lying levels, to allow for cascading 
processes into the corresponding upper levels at n = 7 (which 
are our present uppermost ones). 

2.4. Understanding the Niv /14058/.16380 line formation 

In the following, we discuss the most important mechanisms 
that explain the presence of emission at Niv/14058, in particu- 
lar the decisive role of mass-loss. Our analysis is based on the 
model-grid as described in Sect. 14.21 and refers to LM C back- 
ground abundances plus a solar ( Asplund et al.ll2005h nitrogen 
abundance, [N] = 7.780 chosen in order to obtain pronounced 
effects. 

All the important levels and the corresponding transitions in- 
volved in the Niv emission problem are summarized in Fig. [1] 
A comparison with the analogous diagram for Niii (Fig. 1 in 
Paper I) shows a number of similarities but also differences. In 
addition to what has akeady been outlined, the upper level of the 
emission line (3d — > 3p) is fed by only weak dielectronic recom- 
bination, with almost no influence on the population of 3d (con- 
trasted to the N III case), and there is no resonance line connected 
to 3d (which turned out to be crucial for N iii). Instead, the lower 
level of /14058, 3p, is connected with the ground-state. Similar 



roughly 0.9 dex above the LMC baseline abundance, 
[N]baseiine = 6.9, following lHunter et alj (l2007h . 
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Fig. 1. Simplified Grotrian diagram displaying the most impor- 
tant transitions involved in the N iv emission line problem. The 
horizontal line marks the Nv ionization threshold. Niv /14058 
is formed by the transition 3d 'D — > 3p 'P", while the absorp- 
tion line at /16380 originates from the transition 3p 'P*' 3s 'S. 
An efficient drain of 3p is provided by the 'two-electron' transi- 
tions 3p 'P" — > 2p^ 'S, 'D. Cascade processes from 3d' 'F" and 
pumping from 2p 'P" are the major routes to overpopulate the 
3d 'D state. See text. 
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Fig. 2. Equivalent width (positive for absorption) of NiV/14058 
(black) and Niv/16380 (red) as a function of Teff. Solid and dot- 
ted curves refer to low M (model series A') with logg=4.0 and 
3.7, respectively, and dashed curves to supergiant mass-loss rates 
(model series 'E'), with log §=4.0. 



to Niii, on the other hand, there are two strong 'two-electron' 
transitions able to drain 3p, via 3p 'P" — » 2p^ 'S, 'D. 

2.4.1. Basic considerations 

In agreement with the results from lHeap et al.l (12006 '). our sim- 
ulations (see Fig. |2|i show that NiV/14058 turns from weak ab- 
sorption (around Teff ~ 37 kK) into weak emission around 
Teif ~ 42 kK, for models with (very) low mass-loss rate and 
logg = 4.0. As usual, we define equivalent widths to be posi- 
tive for absorption and to be negative for emission lines. We find 
a 2 kK difference w.r.t. the turning point, which can be attributed. 
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Fig. 3. Departure coefficients, fractional net rates and line pro- 
files for Niv/14058//16380 and involved processes, for model 
'A4540'. Designation for 2p^ refers to 2p^ ^D only since the 
other fine-structure component, 2p^ ^S, behaves similarly. Upper 
panel: NLTE departure coefficients as a function of tross- The 
onset of the wind is clearly visible at tross^ 0.003. Black curves: 
standard model with > b^p; red curves: draining transitions 
3p —> 2p^ suppressed, leading to bj^ < b^^; blue curves: transi- 
tion 2p — > 3d suppressed. Middle panel: fractional net rates to 
and from 3p, for the standard model (black) and the model with 
suppressed draining transitions (red). Lower panel: line profiles 
for N iV/14058 and /16380, for the three models displayed in the 
upper panel, with similar color coding. 



to a major part, to the lower nitrogen content of our models and 
different background abundances. Towards hotter temperatures, 
the emission strength increases monotonically until a maximum 
around Tetf « 53 kK has been reached, after which the emission 



stabilizes and finally decreases. For lower gravities and/or higher 
mass-loss rates, the line turns into emission at lower Teff , so that, 
for a given Teff, the emission strength increases with decreas- 
ing log g and increasing mass-loss rate, M. Since the absorption 
strength of N iv/l6380 increases in a similar way (though with a 
much weaker impact of M, and only until Tgff ~ 50 kK), both 
lines appear (for a given T^s) as anti-correlated, at least for a 
large range of temperatures. In contrast, the corresponding tran- 
sitions of N m were found to be correlated (Paper I). 

The behaviour of both lines and the corresponding level 
structure implies an efficient drain of level 3p that enhances the 
emission at /14058 as well as prevents emission/increases absorp- 
tion at /16380, and is provided by the two 'two-electron' transi- 
tions 3p 2p^ 'D, 'S, similar to the case of N niF^ 

To investigate this mechanism in more detail, and to avoid 
'contamination' by wind effects, we concentrate at first on 
a low-M model with decent emission at NiV/14058, with 
r^ff = 45, 000/: and log ^ = 4.oHln the upperpanel of Fig.[3] 
we provide the NLTE departure coefficients, b, of involved lev- 
els, where black curves refer to our standard model. Obviously, 
level 3d is overpopulated with respect to 3p (/73d > fesp) over the 
complete line formation region. On the other hand, levels 2p^ 'D, 
2p^ 'S, and 2p (the latter two not displayed) are (mostly colli- 
sionally) coupled to the ground state, which in itself is strongly 
depopulated, due to rather large ionizing fluxes (see below). This 
situation closely resembles the situation in N m, where strongly 
depopulated draining levels (for non-blocked models) favoured 
a depopulation of the analogous level 3p. 

To further clarify the impact of the different processes, we 
investigate the corresponding net rates responsible for the popu- 
lation and depopulation of level 3p (Fig. [3] middle panel, black 
curves). As in Paper I, we display the dominating individual net 
rates (i.e., «//?/,■ - niRij > for population, with index / the con- 
sidered level) as a fraction of the total population rate. Indeed, 
the drain by level 2p^ 'D (dashed-dotted) and/or level 2p^ 'S 
(not displayed) are the most important processes that depopulate 
level 3p in the line formation region. In contrast, the resonance 
line does not contribute to any (de-)population of level 3p, since 
it is (almost) in detailed balance (long dashed line). 

To check the validity of our scenario, we calculated an alter- 
native model where the two draining transitions have been sup- 
pressed, by using very low oscillator strengths. Indeed, the upper 
panel of Fig. [3] (red curves) shows that now b^ is smaller than 
fe3p, and /14058 goes into absorption (lower panel, red color). 
From the fractional net rates, we see that the preferred decay 
route has switched from 3p — » 2p^ (standard model, black) to 
3p — > 3s (red, dashed), though level 3p retains a much larger 
population. 

The upper level of /14058, 3d, is predominantly fed by cas- 
cading from 3d' 'F°, and also by pumping from level 2p, whilst 
dielectronic recombinations are negligible. Suppressing the pop- 
ulation from level 2p leads to less emission (Fig.|3] lower panel, 
blue colors), due to a less populated level 3d (upper panel). 

Let us now consider the behaviour of the absorption line at 
/16380, resulting from the transition 3p — > 3s, again by means of 



' ' though in Paper I we argued that in case of N iii this mechanism be- 
comes suppressed in realistic model atmospheres with near-solar back- 
ground abundances. 

At cooler T^g, this line is in absorption because of a lower Niv 
ionization fraction implying deeper formation depths, which are closer 
to LTE. 
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Fig.|3] As mentioned earlier, this line shows an anti-correlatior(3 
with NiV/14058, in contrast to the behaviour of the correspond- 
ing N m lines which appear as correlated. In Paper I we argued 
that the latter correlation results from the proportionality of the 
level populations of 3p and 3s. That is, when b^p decreases (e.g., 
due to increased 'two-electron' drain), decreases in paral- 
lel due to less cascading, and the absorption at /14097 becomes 
weaker in concert with an increase in the triplet emission. Vice 
versa, an increase of 3p implies less emission of the triplet lines 
and more absorption at A4Q91, respectively. 

Such a reaction requires the transition 3p — > 3s to be opti- 
cally thin, dominated by spontaneous decays, which is no longer 
true for NiV/16380. Due to a mostly significant optical depth, 
the radiative net rate is no longer dominated by spontaneous de- 
cays, but depends also on absorption and induced emission pro- 
cesses. Now, an increased population of 3p leads to less increase 
of /73s0 and the absorption becomes weaker because of an in- 
creased source function oc b^p/b^s- Vice versa, a decrease in the 
population of 3p leads to more absorption at /16380 in parallel 
with more emission at /14058. This behaviour becomes particu- 
larly obvious if we investigate the reaction of the absorption line 
when suppressing the draining transitions. In this case, 3p be- 
comes strongly overpopulated (Fig. |3] upper panel, red color), 
and /14058 goes into absorption whilst /16380 becomes an emis- 
sion line, due to a significantly increased source function (more 
pumping than in the original scenario). We checked that if the 
absorption and stimulated emission terms in the 3p —> 3s transi- 
tions are neglected, /16380 displays more absorption instead, in 
accordance with our previous arguments. 

Note, however, that this anti-correlation is not complete. 
E.g., if one changes processes which have an effect on 3d alone, 
the absorption strength of /16380 remains unaltered. Thus, by 
suppressing 2p —> 3d, only b^d is affected (upper panel, blue vs. 
black curves), and there is less emission at /14058 while the ab- 
sorption at /16380 remains at the previous level (lower panel). 

Summarizing, we interpret the different correlations between 
emission and absorption line-strength in N iii and N iv as due to 
optical depth effects in the 3p — > 3s transition. As long as this is 
optically thin, cascade effects dominate, and both lines appear as 
correlated (Nm). Larger optical depths introduce a counteract- 
ing 'source-function effect', and the lines become anti-correlated 
(Niv). The rather large degree of such anti-correlation supports 
the importance of the draining transitions, since these are able 
to influence both the absolute population of the involved levels 
as well as their ratios in a very efficient way, by providing addi- 
tional decay channels for level 3p. 

2.4.2. The impact of wind effects 

So far, we discussed the possibility of obtaining emission at 
NiV/14058 via solely photospheric NLTE processes. In Paper I, 
the presence of a wind (actually, a steep rise of the velocity 
field in the outer photosphere) turned out as crucial to explain 
the observed Nm triplet emission in Of-stars, enabling an ef- 
ficient pumping of the upper level by the corresponding reso- 
nance line. To investigate how the presence of a wind affects the 



when certain parameters/processes are changed /or a given T^ff, 
e.g., the strength of the draining levels, the background opacities and so 
on. The overall increase of these line-strengths as a function of Teff is 
related to the increasing ionization fractions. 

The net radiative rate (downward) is proportional to A(l - J/S), 
with Einstein coefficient for spontaneous decay, A, scattering integral / 
and source function, S . 
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Fig. 4. Departure coefficients (top) and line profiles for 
Niv^4058/ ^6380 (bottom) for the dense-wind model 'E4540'. 
Level designations, lines, and color coding as in Fig.|3] 

emission at N iv/l4058, we compare our previous model 'A4540' 
with model 'E4540', which has the same stellar parameters but 
a considerably larger, supergiant-like mass-loss rate. It turns out 
that the inclusion of such a strong wind has a pronounced ef- 
fect. Comparing Fig. |4](lower panel) with Fig. [3] model 'E4540' 
(black) results in much more emission than 'A4540', increasing 
the equivalent width of /14058 from -7 to -114 mA. Also for 
higher Tgff the impact of M remains significant. The absorption 
line NiV/16380 is affected by the wind as well, though less pro- 
nounced. The absorption becomes slightly stronger (by 18% in 
the equivalent width), i.e., the anti-correlation discussed above 
is still present. This is valid not only for model 'E4540', but also 
for hotter models (Fig.|2ll, until the wind-emission begins to con- 
taminate N iv/i6380. 

The origin of such stronger emission at /14058 becomes 
clear if one inspects the involved departure coefficients (Fig. |4] 
upper panel, black curves). Again, the onset of the wind is 
clearly visible (at tross ~ 0.1), now much deeper than in the 
'A' model, and the line formation region is located in between 
tross ~ 0.40. . .0.04. Compared to the 'A model, 'E4540' dis- 
plays a larger ground-state depopulation, where the ground-state 
remains coupled with the draining levels 2p^ as well as with level 
2p (not displayed). This leads, particularly in the transition re- 
gion between photosphere and wind, to an extreme depopulation 
of 3p. Since level 3d becomes strongly overpopulated, mostly 
due to feeding by 3d' 'F*' (which is severely overpopulated as 
well), the resulting line source function is quite large and partly 
even in inversion, which explains the pronounced emission at 
NiV/14058. 
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All these differences are caused by the onset of the wind. At 
first note that the N iv continuum (/I < 160 A) is strongly coupled 
with t he He ii continuum. As already realized by Gableretal] 
(Il989l) . increasing mass-loss leads to more He n ground-state de- 
population, to higher fluxes in the He ii continuum and thus also 
to higher fluxes in the N iv continuum. Consequently, the N iv 
ground-state becomes strongly depopulated, and the nitrogen 
ionization equilibrium switches from N iv (being the dominant 
stage in the deeper photosphere) towards N v from the transition 
region on. This also favours the overpopulation of 3d, by means 
of increased recombinations to high-lying levels with subsequent 
cascades via 3d' 'F". 

One might now argue that the inclusion of the wind could 
amplify the impact of the resonance transition(s), by produc- 
ing deviations from detailed balance leading to strong pumping, 
similar to the case of N m. Unlike the situation in N m, how- 
ever, such an effect would lead to less line emission or even 
absorption at NiV/14058, since here the resonance line is con- 
nected to the lower level, 3p. As it turns out, however, this effect 
is not present, since the resonance line is too strong (N iv is the 
dominant ion until the transition region) to leave detailed bal- 
ance before the wind has reached a significant speed. Only then, 
the resonance transition becomes dominant in populating level 
3p, but this occurs already far beyond the formation region of 
NiV/14058. Close to the formation region, there is only a mod- 
erate population of 3p by the resonance line, similar to the pop- 
ulation from 3d itself. Even this additional population is coun- 
teracted (actually, even slightly overcompensated) by enhanced 
drain towards 2p^, not only in the formation region of /i4058 
but also in those outer regions where the resonance line strongly 
pumps. 

These arguments are supported by the fact that N iv/i6380 is 
only slightly affected by the wind, where the increased absorp- 
tion results mostly from a diminished source function due to a 
less populated 3p level. Performing the same tests as for the thin 
wind case, i.e., either suppressing the drain or suppressing the 
population of 3d via 2p — » 3d leads to similar results, as can be 
seen from the red and blue curves and profiles in Fig. H) respec- 
tively. 

We conclude that M is a key parameter for modeling the N iv 
emission line, where in contrast to N iii the basic mechanism (for 
typical mass-loss rates and below) is always due to the depop- 
ulation of the lower level by the 'two-electron' transitions. This 
drain becomes stronger as a function of M, because of increasing 
ionizing fluxes leading to more ground-state depopulation. 

3. Stellar sample and observations 

3.1. The stellar sample 

Table [3] lists our stellar sample which has been drawn from the 
analysis of LMC 0-/early B-stars by Mok07. 

Three of the 28 stars from the original sample have been dis- 
carded from the present analysis, for two reasons. First, from 
our analysis we suspect that Nl 1-004 and Nl 1-048 might be 
(SBl) binaries, where the former object shows discrepant line 
shifts and for the latter we were not able to reproduce accu- 
rately the observed He lines (shape and strength). A possible 
binarity of Nil -048 was also suggested by Mok07, because of 
similar reasons. The other discarded star, Sk-67° 166, is the 
only object in the original sample which seems to be strongly 
evolved (helium con tent ^He = N(He)/N(H) = 0.20 . . . 0.28, 
ICrowther et al.ll2002l and Mok07, respectively), and has a very 
dense wind, with both H„ and He ii4686 in strong emission. We 



confirm the stellar/wind parameters as derived by Mok07 (al- 
most perfect fit quality of H/He lines), but did not succeed in 
a reasonabl e fit for the nit r ogen lines. A comparison with the 
analysis bv ICrowther et al.l (l2002l) shows similar discrepancies. 
Because of this problem and because of its highly evolved evo- 
lutionary status, which does not match with all other objects in 
our sample, we decided to discard this object from our present 
analysis and will re-consider it in a future attempt. 

The remaining sample consists of 20 O-stars, mostly giants 
or dwarfs, and 5 early B-type supergiants or giants. All of the B- 
stars and 15 O-stars are associated with the cluster Nil, and the 
others are field stars. The early B-stars have been included in our 
sample to allow us for a com parison with previous analyses of 
such stars dHunter et al.ll2009h in order to check the consistency 
of different codes and methods in the transition region between 
O- and B-types (Sect.lOli. 

Table[3]gives information about spectral type, V-magnitude, 
interstellar extinction and absolute visual magnitude, and has 
been tak en from Mok07. S pectral types for Nil objects are 
based on lEvans et all (|2006|) . slightly revised by Mok07 in col- 
laboration with C. Evans (priv. comm.), and for the field stars 
from Walbor net af] (Il995h . and lMassev eTaP (1 19951 l2005h . For 
the field star Sk-70° 69 we added the ((f)) designation since the 
present spectra show clear emission at N iuAA4634- - 4640-4642 
and He ii/t4686 in absorption. 

3.2. Observations 

Most of the observations (for objects denoted by 'Nil-') have 
been carried out within the FLAM ES I survey, and are de- 
scribed in detail in lEvans et al.l (l20d^. In brief, the data were ob- 
tained using the Fibre Large Array Multi-Element Spectrograph 
(FLAMES) at the VLT, for six wavelengths settings with an ef- 
fective resolving power of R ^ 20, 000. The S/N ratios are 
in the ranges 50-200 for LMC objects. After sky substraction, 
each wavelength range was co-added and normalized by means 
of a cubic spline[3 The final merged spectra cover two spectral 
ranges, 3850-4750 A and 6300-6700 A. 

To improve the sampling in luminosity and temperature, 
Mok07 augmented the Nil sample by LMC O-type field stars, 
which were observed using the UVES spectrograph at the VLT 
as part of the ESO programmes 67.D-0238, 70.D-0164, and 
074.D-0109 (P.I. P. Crowther). Spectra were obtained for four 
different wavelength settings, at an effective resolving power of 
R - 40, 000. The final product provides coverage between 3300- 
5600 A and 6300-10400 A for all stars except for Sk-70° 69, 
where 'only' the region between 3300-5600 A and 6300-6700 A 
had been observed. The typical S/N ratios achieved for all spec- 
tra lie in the range 60-80. 



4. Analysis 

4.1. Methodology 

In an ideal world, we could have used the stellar and wind pa- 
rameters as obtained by Mok07 from H/He lines, and simply 
derived the nitrogen abundances on top of atmospheric models 
with these parameters. Unfortunately, there are reasons to re- 
analyze all programme stars. At first, the present fastwind ver- 
sion (see Sect. 12.1b is somewhat different from the version used 



We performed additional re-normalizations for different wave- 
length ranges. 
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Table 3. Sample stars used within this study, along with spectral 
type, V-magnitude, interstellar extinction A v, and absolute visual 
magnitude My- All quantities have been taken from Mok07. 



Star 


Cross-IDs 


Spectral Type 


V 


Av 


Mv 


Nil -026 


_ 


02 111(f) 


13.51 


0.47 


-5.46 


Nll-031 


P3061/LH10-3061 


ON2 111(f) 


13.68 


0.96 


-5.78 


Nll-038 


P3100 


05 ll(f+) 


13.81 


0.99 


-5.68 


Sk-66° 100 




06 11(f) 


13.26 


0.34 


-5.58 


Nl 1-032 


P3168 


07 11(f) 


13.68 


0.65 


-5.47 


Nl 1-045 




09111 


13.97 


0.50 


-5.03 


B1253 




02 V((f )) 


13.76 


0.71 


-5.45 


B1237 


- 


02 V((f )) 


13.89 


0.62 


-5.23 


Nl 1-060 


P3058/LH10-3058 


03 V((f )) 


14.24 


0.81 


-5.07 


Sk-70° 69 


_ 


05 V((f)) 


13.95 


0.28 


-4.83 


Nl 1-051 


_ 


05 Vn((f)) 


14.03 


0.19 


-4.66 


Nl 1-058 


_ 


05.5 V((f)) 


14.16 


0.28 


-4.62 


Sk-66° 18 


_ 


06 V((f)) 


13.50 


0.37 


-5.37 


Nl 1-065 


PI 027 


06.5 V((f)) 


14.40 


0.25 


-4.35 


Nl 1-066 




07 V((f)) 


14.40 


0.25 


-4.35 


Nl 1-068 




07 V((f)) 


14.55 


0.28 


-4.23 


Nl 1-061 




09 V 


14.24 


0.78 


-5.04 


Nll-123 




09.5 V 


15.29 


0.16 


-3.37 


Nl 1-087 


P3042 


09.5 Vn 


14.76 


0.62 


-4.36 


Nl 1-029 




09.7 lb 


13.63 


0.56 


-5.43 


Nl 1-036 




B0.5 lb 


13.72 


0.40 


-5.18 


Nl 1-008 


Sk-66° 15 


B0.7 la 


12.77 


0.84 


-6.57 


Nl 1-042 


P1017 


BO 111 


13.93 


0.22 


-4.79 


Nl 1-033 


PI 005 


BO llln 


13.68 


0.43 


-5.25 


Nl 1-072 




B0.2 111 


14.61 


0.09 


-3.98 



Note s. Primary identifications for Nil objects are from lEvans et alJ 
( 1200 6). Identifications starting with "Sk", " Bl", "P", and "LH " are 
from Sanduleak ( 1970), Brun et et alJ dl975h . iParker et alj ( Il992h . and 
IWalborn et alj (i2002i.i2004) . respectively. 



by Mok07, and the parameters need certain (mostly small) al- 
terations in order to reach a similar fit quality to the H/He lines. 
Second, we used a somewhat different fitting strategy with re- 
spect to the 'free' parameters, which changes the optimum fit. In 
contrast to Mok07, we derived and fixed the projected rotational 
velocity, v sin ;, independently from the actual fitting procedure 
(Sect. 14.3b . whereas Mok07 included v sin / as a free parame- 
ter in their fitting algorithm. Moreover, during our fit procedure 
we allowed for the presence of extra line-broadening ('macro- 
turbulence', Vmac). not considered by Mok07. Differences in 
V sin / and Vmac can lead to certain differences in the outcome 
of the fit, since the profile shapes might change (e.g.. Fig. 4 in 
lPulsll2008 t). Third, and most important, is the fact that we now 
aim at a consistent fit for the H/He and nitrogen lines. Thus, and 
in the sense of a compromise solution (minimization of the dif- 
ferences between observed and synthetic spectra for all lines), 
different stellar parameters which result in a modest change of 
the fit quality of H/He alone[]3 can lead to a significant improve- 
ment with respect to the complete set of lines. 

In so far, we opted for an entire re-analysis, performed 
mostly by a simple 'fit-by-eye' method where we try to accom- 
plish the best fit to the strategic lines by visual inspection (for 
a discussion, see i Mokiem et alj|2005h . Since we start from the 
parameter set as provided by Mok07 (highest 'fitnessQ with re- 

In particular at earliest spectral types, the sensitivity of H/He on 
changes in the atmospheric parameters is rather weak. 

" which quantifies the quality of the solutions resul ting from a ge- 
netic algorithm optimization, see iMokiem et al.l ( 120051) and references 
therein. 



spect to their assumptions), our new solution should be located 
at or close to the global maximum of the corresponding merit 
function as well and not only at a local one. Note that our deriva- 
tion of nitrogen abundance and micro-turbulence for most of the 
cooler sample stars relies on a more objective method (Sect. l4.4T l. 

MokO ? themselves used an automated fitting method (devel- 
oped by Moki em et al.l l2005) based on a genetic algorithm opti- 
mization routine to obtain the stellar/wind parameters by evolv- 
ing a population of fastwind models over a course of genera- 
tions, until the best fit to H/He is found. Seven free parameters 
were considered to obtain the highest 'fitness': effective temper- 
ature, Tf.fl, surface gravity, \ogg, helium content. Fee. projected 
rotational velocity, v sin / (see above), micro-turbulent velocity, 
Vmic. mass-loss rate, M, and velocity-field exponent,/?. 

4.2. Model calculations and grids 

All models used within this analysis are calculated with fast- 
wind, augmented by few cmfgen models for comparison pur- 
poses (Sect.|5]l. For these calculations, H, He and N are treated 
as explicit elements. A description of our H/He model atoms can 
be found in Puis et al. (2005), and our nitrogen model atom has 
been described in Sect. l2.2l and in Paper I. 

To allow us for studying the combined reaction of all di- 
agnostic H/He/N lines on variations of the stellar/wind param- 
eters and nitrogen abundances, and also for understanding the 
NiV/}4058 emission line process (Sect l2.4l ). we generated a grid 
of modelsHThe grid is constructed using various nitrogen abun- 
dances centered at the solar value [N] - 7.78 (from [N] - 6.98 
to 8.58 with step size 0.2 dex and including the LMC nitrogen 
baseline abundance, [N]baseiine - 6.9), and a background metal- 
licity, Z = 0.5Zq, correspoi iding roughly to the g lobal metallic 
abundance of the LMC (cf. lMokiem et al.ll2007bl) . The individ- 
ual abundances of the background elements (in terms of mass 
fractions) are scaled by the same factor with respect to the solar 
abundance pattern (see iMassev et al.ll20041) . 

For a given Z, the grid is three-dimensional with respect to 
Teff , log g, and log Q, where Q is the so-called wind-strength pa- 
rameter (or optical depth invariant), Q = M/(voo/^.)'^. This 
parameter allows us to condense the dependence on M, termi- 
nal velocity, Vco, and stellar radius, R,, into one representative 
quantity^ The grids cover the temperature range from 25 to 
55 kK (with increments of 1 kK), and a gravity range between 
3.0 and 4.5 (with increments of 0.2 dex). For log Q, the differ- 
ent wind strengths are denoted by a letter (from 'A' to 'E'), with 
logg = -14.0, -13.5, -13.15, -12.8, -12.45, respectively, if 
M is calculated in Moyr"', Voo in kms"', and Rt in Rq. Models 
with quantifier 'A' correspond to thin winds, resulting in lines 
that are (almost) unaffected by the wind, whereas 'E' -models 
correspond to a significant wind-strength typical for O-type su- 
pergiants. Other parameters have been adopted as follows: a so- 
lar helium abundance, iee = 0.10; Vm as a function of the pho- 
tospheric escape velocity, Vesc (see Kudritzki & Puis 2000); the 
stellar radius, Rt, as a function of spectral type and luminosity 
class, corresponding to prototypical values; the velocity field ex- 
ponent, /3, from empirical values (Kudritzki & Puis 2000), with 
/3 = 0.9 for O-stars, and higher values towards later types; and 
the micro-turbulence, Vj^ic = 10 kms"'. 



'** following the basic philosophy described bv lPuls et al .' ('2005). 

" relying on the fact that the wind-emission from recombination dom- 
inated (i.e., p--) processes remains unaffected as long a s the wind- 
strength parameter does not vary fsee lPuIs et al.ll 199^120051) . 
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Table 4. Fundamental parameters for the LMC sample, assuming unclumped mass-loss. 



Star 


ST 


Secondary 






log^ftiue 


R. 


log L. 






V sin i 




M 


P 


I'co 






T^ff diag. 


(kK) 


(cgs) 


(cgs) 


{Re) 


(Le) 




(km s"') 


(km s-') 


(km s" 


') (lO-^Moyr- 




(km s-') 


Nl 1-026 


02 m(f ) 


Niii/Niv/Nv 


49.0 


4.00 


4.00 


11.3 


5.82 


0.10 


10.0 


72 


60 


1.56 


1.08 


[3120] 






Niv/Nv 


52.0 


4.10 


4.10 


11.0 


5.89 










1.49 






Nll-031 


ON2 111(f) 


N m/N IV 


47.8 


3.95 


3.95 


13.4 


5.92 


0.11 


10.0 


71 


60 


2.02 


1.08 


3200 






Niv/Nv 


56.0 


4.00 


4.00 


12.2 


6.12 










2.20 






Nl 1-038 


05 ll(f+) 


N iii/N IV 


40.5 


3.70 


3.71 


14.0 


5.67 


0.08 


10.0 


100 


50 


1.21 


0.98 


[2600] 


Sk-66° 100 


06 11(f) 


N m/N IV 


39.0 


3.70 


3.71 


13.7 


5.59 


0.19 


10.0 


59 


60 


0.83 


1.27 


2075 


Nl 1-032 


07 11(f) 


N m/N IV 


36.0 


3.50 


3.51 


13.6 


5.44 


0.09 


10.0 


60 


70 


0.97 


0.80 


[1920] 


Nl 1-045 


09 111 


- 


32.3 


3.32 


3.33 


12.0 


5.15 


0.07 


10.0 


64 


80 


0.69 


0.80 


[1550] 


B1253 


02 V((f )) 


Niv/Nv 


54.8 


4.18 


4.20 


10.7 


5.97 


0.08 


10.0 


230 


- 


1.53 


1.21 


3180 


B1237 


02 V((f )) 


Niv/Nv 


53.2 


4.11 


4.12 


9.7 


5.83 


0.09 


10.0 


140 


- 


0.62 


1.26 


3400 


Nl 1-060 


03 V((f )) 


Nm/Niv/Nv 


48.0 


3.97 


3.97 


9.5 


5.63 


0.12 


10.0 


68 


40 


0.51 


1.26 


[2740] 






Niv/Nv 


51.0 


4.10 


4.10 


9.2 


5.71 










0.48 






Sk-70° 69 


05 V((f)) 


N m/N IV 


42.3 


3.93 


3.94 


9.1 


5.38 


0.14 


10.0 


131 


- 


0.43 


0.80 


2750 


Nll-051 


05 Vn((f)) 


N m/N IV 


41.4 


3.70 


3.83 


8.6 


5.42 


0.08 


10.0 


350 


- 


0.41 


0.80 


[2110] 


Nl 1-058 


05.5 V((f)) 


N m/N IV 


40.8 


3.75 


3.76 


8.4 


5.24 


0.10 


6.0* 


62 


60 


0.01 


1.00 


[2470] 


Sk-66° 18 


06 V((f)) 


N m/N IV 


39.7 


3.76 


3.76 


12.2 


5.52 


0.14 


10.0 


75 


40 


1.08 


0.94 


2200 


Nl 1-065 


06.5 V((f)) 


N m/N IV 


41.0 


3.85 


3.85 


7.4 


5.14 


0.13 


10.0 


60 


50 


0.05 


1.00 


[2320] 


Nl 1-066 


07 V((f)) 


N m/N IV 


37.0 


3.70 


3.71 


7.9 


5.02 


0.10 


10.0 


59 


50 


0.14 


0.80 


[2315] 


Nl 1-068 


07 V((f)) 


N m/N IV 


37.0 


3.70 


3.71 


7.5 


4.98 


0.10 


10.0 


30 


50 


0.13 


1.12 


[3030] 


Nl 1-061 


09 V 


- 


34.0 


3.55 


3.55 


11.5 


5.20 


0.09 


5.0 


54 


80 


0.52 


0.80 


[1900] 


Nll-123 


09.5 V 




34.3 


4.20 


4.21 


5.4 


4.56 


0.09 


10.0 


115 




0.08 


0.80 


[2890] 


Nl 1-087 


09.5 Vn 




32.7 


4.04 


4.08 


8.8 


4.90 


0.10 


10.0 


260 




0.11 


0.80 


[3030] 


N 11-029 


09.7 lb 


Nii/Niii 


29.0 


3.20 


3.21 


15.8 


5.20 


0.08 


10.0 


46 


60 


0.28 


1.23 


[1580] 


Nl 1-036 


B0.5 lb 


Nii/Niii 


25.8 


3.11 


3.11 


15.5 


4.98 


0.08 


12.1* 


39 


40 


0.11 


0.80 


[1710] 


Nl 1-008 


B0.7 la 


Nii/Niii 


26.3 


3.00 


3.00 


29.5 


5.57 


0.10 


4.7* 


46 


60 


0.62 


1.30 


[2390] 


Nl 1-042 


BO 111 


Nn/Nm 


29.2 


3.59 


3.59 


11.9 


4.97 


0.08 


5.0 


21 


25 


0.19 


1.19 


[2310] 


Nl 1-033 


BO llln 




26.7 


3.20 


3.34 


15.7 


5.05 


0.10 


5.0 


256 




0.25 


1.03 


[1540] 


Nl 1-072 


B0.2 111 


N ii/N III 


29.8 


3.70 


3.70 


8.1 


4.67 


0.10 


2.7* 


14 


10 


0.25 


1.30 


[2100] 



Notes, log^tme is the surface gravity corrected for centrifugal effects. Secondary T^ff diagnostics, used in parallel with the Hei/Heii ionization 
equilibrium, is indicated for each star. For three stars, two parameter sets are provided, due to problems of reaching a simultaneous fit for all 
considered lines from three ionization stages. For Nl 1-026 and Nl 1-060, the first entry is our preferred one, and the second entry is a (disfavoured) 
alternative, whilst for Nll-031 we consider both solutions as possible, (see Sect.|5j. Brackets around v„: derived from scaling via v^^c- Asterisks 
following Vmic! value obtained in parallel with nitrogen abundance. For errors, see text. 



4.3. Determination of stellar and wind parameters 

The different steps performed in our analysis can be summarized 
as follows, and are detailed in the next sections. 

First, we determine v sin / for each object. Then we use 
the results from Mok07 in combination with our model grid to 
roughly constrain the stellar (Tss, logg, [N]) and wind-strength 
parameters (log Q), by inspecting the synthetic and observed 
H/He and nitrogen line profiles (Yue, /3 and v^ic already specified 
within the grid.) During this step, we determine the extra line- 
broadening parameter, v^ac, by reproducing the profile shape of 
the weaker lines. 

Subsequently, the stellar/wind parameters, now including 
Fhc and /3, are fine-tuned by calculating a grid of much higher 
resolution around the initial guess and adopting Voo from Mok07. 
After the fundamental parameters have been fixed and, in case, 
Vmac has been adjusted, we fit the nitrogen abundance, using two 
different methods (for cooler and hotter objects, respectively), 
and also update Vmic, which Mok07 solely derived from H/He 
lines. In certain cases, we need to re-adapt the stellar/wind pa- 
rameters to obtain the (almost) final solution. Now, we are able 
to calculate the stellar radius from My and the synthetic fluxes, 
and we can update the mass-loss rate to its final value, by scaling 
with the new radius. A final consistency check with the new M 
and values is performed to ensure the stability of our results. 

Tables |4] and |6] list all quantities derived in this way, and 
Table |5] yields the main differences between our and the Mok07 



results. Note that these quantities refer to unclumped winds, 
whilst in Sect. l4.5.T] we discuss the impact of wind clumping. 

Projected rotational velocities and macro-turbulence. Before 
we are able to perform the actual (fine-)analysis, we need to 
constrain the line broadening parameters, i.e., v sin / and Vmac 
(v'mic will be - when possible - inferred in parallel with the ni- 
trogen abundances, see Sect l4.4l ). As outlined above, Mok07 de- 
rived V sin / directly from their automated fitting method, from 
the H/He lines. More suitable is to use metal lines, since these 
are not affected by Stark-broadening. Only in case of high ro- 
tational velocities or high temperatures, where metallic lines 
are blended or are very weak. He i lines might be used. Thus, 
we derived v sin / from scratch, employing the Fourier method 
(Grav 1976), as implemented a nd tested in the OB-star rang e 
by Simon -Diaz et al.l (12006) and Simon-Dfaz & Herrerol (120071) . 
This method has the advantage to easily discriminate the rota- 
tional contribution from other broadening mechanisms that af- 
fect the line shapes. In dependence of temperature and rotational 
velocity of the star, we used lines from O ii, N ii, C ii, and Si iii 
for B- and late O-type stars. For earlier O-types, higher ioniza- 
tion states are predominant, and mostly N iii/N iv and Si iv lines 
have been considered, together with He i lines for the earliest 
types. 

To finally reproduce the actual profile shape, in most cases 
some extra line-broadening is needed, conventionally called 
macro-turbulence. Though the physical origin of this broaden- 



11 



J.G. Rivero Gonzalez et al.: Surface nitrogen abundances for 0-stars in the Large Magellanic Cloud 



ing remains still to be proven, there are some strong indications 
that it is associated with (high order, non-radial) stellar pulsa- 
tions ([Xerts et al. 2009; Simon-Diaz et al. 2010). To account for 
this effect, we used a radial-tangential description of Vmac to fit 
the profile shapes of nitrogen (and partly helium) lines, using 
our first estimates on the stellar and wind parameters (see above) 
and the new v sin / values. For the fastest rotators of our sample, 
however, corresponding values could not been constrained, since 
large v sin / produce either too weak nitrogen lines, or these 
lines, together with He i lines, loose their sensitivity to distinct 
changes in v'mac- 

As expected, most of our v sin / values turn out to be system- 
atically lower than those provided by Mok07 (Table |5]l, by typi- 
cally 30-40%. The derived range of v„iac va lues is consist e nt wit h 
results from similar investigations, e.g., Dufton et al.' (2006*); 
[Simon-Diaz et al. (2006); Lefever et al. (2007); Markova & Puis 
( l2008b ; ISimon-Diaz et all (12010) . The uncertainty of our esti- 
mates is typically on the order of ±10 kms ^ being larger for 
stars with relatively low rotational speeds. 



Effective temperatures. This parameter is mostly constrained 
by the He i/He ii ionization equilibrium. For this purpose, we pri- 
marily use UeiAA 4471, 4713, 4387 and Hen AA 4200, 4541. In 
most cases, we did not meet the so-called He i singlet problem 
dNaiarro et alJl2006l) . As a consistency check on T^fi and espe- 
cially for the hotter stars, where the He i lines can no longer serve 
as an efficient temperature indicator, we make additional use of 
the nitrogen ionization equilibrium, by means of the lines listed 
in Table |2] For B- and late O-stars, we investigate N ii/N iii, for 
mid O-stars N iii/N iv, and for early O-stars N iv/N v or even - 
in a few cases -N iii/N iv/N v (see Table |4] for the specific diag- 
nostics applied to a particular object). To this end, we either use 
our coarse grid or our specific 'fine-grid' models, with a sim- 
ilar gridding of nitrogen abundances ([N] = 6.9 . . . 8.58). By 
exploiting this additional information, i.e., roughly 'fitting' the 
nitrogen lines of different ions at a unique abundance, we are 
able to fine-tune T^ff (and also some of the other parameters, see 
below). For objects where only lines from one ionization stage 
are present (Nl 1-033, 045, 061, 087, 123), such a concistency 
check only allows for rather weak constraints, if at all. 

We estimate a typical uncertainty for T^ff according to the 
grid resolution, A^eff ~ 1 kK. For Nl 1-066 and Nl 1-068 we 
can only provide rough estimates on the stellar parameters, con- 
sistent with the nitrogen ionization equilibrium, and we adopt a 
larger error, AFeff ^ 2 kK. For some problematic stars, Nl 1-026, 
Nl 1-031, and Nl 1-060, we consider an even larger eiTor, about 
4 kK, in agreement with Mok07 (see Sect.|5]l. 

Surface gravities. We derived \ogg using the classical ap- 
proach, from the Stark-broadened wings of the Balmer lines, 
basically Hy and H^, which should be uncontaminated by wind- 
emission. As for Tes, we use nitrogen lines as a final consistency 
check. These surface gravities need to be coiTected for stellar ro- 
tation, applying a centrifugal correction (see Repolust et al. 2004 
and references therein). The estimated error for \ogg is 0.1 dex. 

Helium abundances. To ensure the reliability of the final pa- 
rameters, especially Tgg and and log g, and for our discussion 
on the abundance enrichment, we need to revisit the helium line 
fits, since inconsistent helium abundances can influence these 
parameters. A small sub-grid was constructed around the stellar 
parameters derived in previous steps, for different Yue, from 0.08 



Table 5. Differences between the fundamental parameters as de- 
rived here and by Mok07. 



Star 




Alogg 


AKho 


Av sin i 


AlogM 


A/3 




(kK) 


(cgs) 




(kms-') 


(10-^Moyr-') 




Nil -026 


-4.3 


- 


-0.01 


-37 


-0.06 


- 


Nll-031 


2.8 


0.10 


0.01 


-45 


-0.28 


0.19 


Nll-038 


-0.5 


-0.02 


-0.02 


-45 


-0.10 


- 


Sk-66° 100 


- 


- 


- 


-25 


-0.03 


- 


Nll-032 


0.8 


0.05 


- 


-36 


0.08 


-0.23 


Nil -045 


- 


- 


- 


-41 


0.10 


- 


BI253 


1.0 


- 


-0.01 


39 


-0.10 


- 


BI237 


- 


- 


-0.01 


14 


-0.10 


- 


Nil -060 


2.3 


0.05 


- 


-38 


-0.01 


- 


Sk-70° 69 


-0.9 


0.06 


-0.03 


- 


-0.38 


0.02 


Nl 1-051 


-1.0 


-0.05 


- 


17 


-0.39 


0.20 


Nll-058 


-0.5 


-0.14 


- 


-23 


-1.18 


-0.42 


Sk-66° 18 


-0.5 


- 


- 


-7 


- 


- 


Nil -065 


-0.7 


-0.04 


-0.04 


-23 


-0.86 


0.20 


Nil -066 


-2.3 


-0.17 


-0.01 


-12 


-0.47 


- 


Nil -068 


-2.9 


-0.43 


- 


-24 


-0.42 


- 


Nll-061 


0.4 


0.04 




-33 


0.39 


-1.00 


Nll-123 


-0.5 


-0.02 




5 






Nil -087 








-16 


-0.10 




Nil -029 


-0.4 


-0.03 


0.01 


-31 


0.22 


-0.40 


Nll-036 


-0.5 


-0.20 




-15 






Nil -008 


0.3 


0.02 




-35 


0.09 


-0.57 


Nil -042 


-1.0 


-0.10 


-0.02 


-21 


-0.02 




Nll-033 


-0.5 


-0.01 


0.02 




0.02 




Nil -072 


-1.0 


-0.08 


-0.02 




0.02 


0.46 



Notes. Positive values indicate larger values from this study. Dashes: no 
difference compared to the Mok07 analysis. For Nll-031, the compar- 
ison is made with respect to the cooler solution (see Sect.O. 



(corresponding to the approximate LMC baseline abundance, 
see Sect. |6]l to 0.14, in steps of 0.02. A rough estimate on the 
error is half this stepsize. 



Terminal velocities cannot be reliably derived from the optical, 
and have been adopted from Mok07. For the field stars, values 
have been inferred from UV P Cygni profiles, by IMassa et al.l 
(2003) and Massey et al. (2005). For Sk-66° 18, it was measured 
by Mok07 using UV O vi lines. For the FLAMES Nl 1 stars, only 
Nll-031 could be an alyzed with respect to this parameter, by 
IWalborn et al ] (I2004|K For all other stars, Vpo has been estimated 
from Vesc, foUowing lKudritzki & PuTsI (|2000|) . 



Velocity field exponent The sample used within this study 
does not contain any star with such a dense wind that Hq, is in 
emission, therefore an accurate determination of this parameter 
is difficulj^ for optical spectroscopy. We applied the following 
philosophy. If the combination M-fi provided by Mok07 resulted 
in reasonable H^-fits, we kept Otherwise, we set /3 to prototyp- 
ical values, /3 = 0.8 ... 1 .30, in dependence of spectral type and 
results from earlier analyses performed in our group. Moreover, 
for some of the sample stars, Nl 1-008, N-029, Nll-058, and 
Nl 1-061, the automated fitting method used by Mok07 resulted 
in quite large values foryS (e.g., N 1 1 -06 1 : y8= 1 . 8), whilst for N 1 1 - 
05 1 a rather low value, /3 = 0.6, was inferred. We consider such 
values as either unphysical or indicating a substantial amount 



^° if not impossible, due to the M- p degen eracy (e.^.. [ Markova et al.l 
l2004h . and effects from wind clumping fe.g.. lPuls et alj|2006h . 
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of wind clumping. In all these cases, we modified j8 as outlined 
above (see also Sect.|5]l. 

Mass-loss rates are derived from fitting the synthetic H^^ pro- 
files to the observations, given (see above). Usually He n/i4686 
is used as a consistency check. It turned out that for many sam- 
ple stars we were not able to successfully fit both lines at the 
same M, since He ii/l4686 showed more absorption than consis- 
tent with the observations when H„ was fitting. This might indi- 
cate a certain problem regarding He ii/l4686 in the new fastwind 
version, or some impact of wind clumping. The problem needs 
to be investigated in future, but has no impact within the present 
study. 

Another consistency check for M is provided by the nitro- 
gen lines, particularly by the N m and N iv emission lines (and 
sometimes also by the N v doublet), which are strongly affected 
by the wind strength. In the case of two stars. Nil -058 and Nl 1- 
065, which showed quite good line fits to H, He and N, the syn- 
thetic N iV/14058 profile displayed weak emission, whereas the 
observed one was clearly in weak absorption. Consistency could 
be achieved by lowering M until this line could be acceptably 
fitted, leaving the remaining nitrogen lines and H„ almost un- 
altered. In both cases it turned out that H„ was already almost 
insensitive to reductions in M. Other stars which showed a simi- 
lar problem, Nll-051, Sk-66° 18, and Sk-66° 100, could not be 
'cured' by this approach, since in this case a reduced value of M 
was no longer consistent with (unclumped!, see Sect. l4.5.TT i H„. 

Since in all of our objects Hq, is in absorption, leading to the 
well known M-p degeneracy, we estimate quite a large error on 
M, namely pl us/minus a factor of t wo, which is typical in this 
situation (e.g.. lMarkova et al.ll2004l) . The impact of the error in 
Rt is negligible here, as outlined in the next paragraph. 

Stellar radii. Since the effective temperatures derived within 
this work are different from those of Mok07 (overall, these dif- 
ferences are modest, except for Nl 1-026, Nl 1-031, and Nll- 
068, see Table |5]l, this leads to different theoretical fluxes and 
thus to di f ferent stell ar radii. Simil a r to M ok07, we followed 
iKudritzkil (Il980h and iHerrero et aTl (Il992h . and calculated the 
'new' radii from the theoretical Eddington fluxes and the (de- 
reddened) absolute magnitudes from Table |3] 

Once the radii have been re-determined, M needs to be mod- 
ified as well, to preserve the fit quality of H„, which depends on 
the optical depth invariant Q (see above). Contrasted to the case 
of Galactic objects, where the eiTor of My (due to unknown dis- 
tances) dominates the eiTor budget of M, this plays a secondary 
role in our sample, due to sufficiently well-known distances and 
the rather large error introduced by the M-f5 degeneracy. 

4.4. Nitrogen abundances and micro-turbulences 

After having determined the stellar and wind parameters and 
their uncertainties (some fine-tuning might still be necessary), 
we are now in a position to derive the nitrogen abundances and 
the corresponding micro-turbulent velocities, v,nic. Because the 
latter parameter significantly affects the strength of both He and 
metal lines and thus the implied abundances, it is useful to de- 
termine both quantities in parallel. 

To caiTy out this analysis, we calculate a fine grid of typically 
25 models, by combining different abundances centered at the 
rough estimates derived in Sect. 14.31 with five different values 
forvmic = 0,5,10,15,20 kms-i. 



Table 6. Helium and nitrogen abundances for the LMC sample, 
with stellar parameters from Table ID 



Star 


1' sin i 




[N] 


A[N] 


A[N]^' 


Literature 


Nil -026 


72 


0.10 


7.80' 


+0.40 


0.25-0.30 


- 








7.75- 


+0.40 






Nll-031 


71 


0.11 


7.83' 


±0.15 


- 


8.00±0.18 








8.30^ 


±0.15 






Nl 1-038 


100 


0.08 


7.85 


±0.15 


0.05-0.10 


- 


Sk-66° 100 


59 


0.19 


8.48^ 


±0.15 


- 


- 


Nl 1-032 


60 


0.09 


7.87 


±0.15 


0.05-0.10 


- 


Nl 1-045 


64 


0.07 


6.98 


±0.20 


- 


- 


BI253 


230 


0.08 


7.90 


±0.15 


0.15-0.20 


- 


BI237 


140 


0.09 


7.38 


±0.15 


0.05-0.10 


- 


Nl 1-060 


68 


0.12 


8.20' 


+0.30 


0.05-0.10 


- 








8.15- 


+0.30 






Sk-70° 69 


131 


0.14 


8.053 


±0.15 


0.15-0.20 


- 


Nll-051 


350 


0.08 


7.58^ 


±0.20 


- 


- 


Nl 1-058 


62 


0.10 


8.09 


±0.15 


- 


- 


Sk-66° 18 


75 


0.14 


8.48^ 


±0.15 


- 


- 


Nl 1-065 


60 


0.13 


8.173 


±0.15 


0.05-0.10 


- 


Nl 1-066 


59 


0.10 


8.17 


±0.20 


0.05-0.10 


- 


Nl 1-068 


30 


0.10 


7.85 


±0.20 


0.05-0.10 


- 


Nl 1-061 


54 


0.09 


7.18 


±0.15 






Nll-123 


115 


0.09 


7.00 


±0.15 






Nl 1-087 


260 


0.10 


7.38 


-0.20'* 






Nl 1-029 


46 


0.08 


7.43 


±0.15 




7.10±0.35 


Nl 1-036 


39 


0.08 


7.85 


±0.17 




7.76±0.11 


Nl 1-008 


46 


0.10 


8.08 


±0.11 




7.84±0.11 


Nl 1-042 


21 


0.08 


7.00 


±0.15 




6.92±0.24 


Nl 1-033 


256 


0.10 


7.28 


-0.20^ 






Nl 1-072 


14 


0.10 


7.68 


±0.15 




7.38±0.06 



From N m/N iv, uncertainty from N v. 
From N iv/N v, uncertainty from N in. 
Improvement of N iV/14058 fit with weak clumping. 
Uncertainty of upper limit. 



Notes. [N] is our best-fitting value, A[N] the corresponding uncer- 
tainty, and A[N]'^' the change in [N] (always positive) if clumping with 
/„ = 0.1 is included in the models (Sect. l4.5"Tt . Values displayed in ital- 
ics have been derive d by the 'curv e-of-growth' method. Literat ure val- 
ues are from Walbo m et al ] ( 120041) for Nll-031 and from Hunt er et al.l 
120091) for the B-stars and Nl 1-029. When two entries are provided, 
the first one corresponds to our preferred solution, except for Nl 1-031 
where we consider both solutions as possible. 



When possible (see below), we use a 'curve of growth' 
method based on the equivalent widths of the lines, which has 
been applied in the past years to different se ts of B-stars for 
obtaining various metallic abundances (e.g., Urb aneial l2004t 
Simon-Diaz et al. 2006; Markova&Puls 2008). In brief, this 
method uses synthetic and observed equivalent widths including 
uncertainties from all considered lines, to derive a unique pair of 
abundance and Vmic (incl. errors). Results from such analyses are 
indicated by an asterisk in the Vn,ic-column of Table|4] 

Unfortunately, this procedure was not applicable to the bulk 
of the sample stars because of various reasons, e.g., blending 
and diluted lines due to fast rotation, almost invisible lines due 
to low nitrogen content, and peculiarities in the observed lines 
from some of the stars with highly ionized nitrogen. All these 
problems will be commented on in Sect. |5] and we opted for a 
determination of the nitrogen abundance/micro-turbulence pair 
by means of a visual inspection of the fit quality, using the same 
fine grid as described above. 

We estimate the errors in Vmic as 3-5 kms"', both for the 
equivalent width and the visual method. To provide an impres- 
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sion on the impact of such errors on the derived nitrogen abun- 
dances, we note that a decrease in Vmic by 5 kms"' leads to 
an increase of [N] by 0.05-0.07 dex. For the error associated 
to [N] when derived by 'visual' fitting, we decided to be quite 
conservative. Even though we are able to obtain quite good fits 
for the bulk of the stars implying an uncertainty of 0.1 dex, we 
rather adopt a larger value of 0. 15 dex to roughly account for the 
additional dependence on the stellar and wind parameters. For 
two stars, Nl 1-033 and Nl 1-087, we can only provide an upper 
limit on the abundance, with an estimated uncertainty of 0.20 
dex. This value is also adopted for Nl 1-066, Nl 1-068, and Nl 1- 
045, Nl 1-051. For the first two stars only rough estimates on 
their stellar parameters could be obtained, and for the latter two 
only one nitrogen multiplet could be used. Even larger uncer- 
tainties have been derived for two 'problematic' stars, Nl 1-026 
and Nl 1-060, see Sect.|5] 

Table |6]lists the obtained nitrogen abundances together with 
their estimated errors. When available, coiTesponding literature 
values have been added. 

4.5. Additional considerations 
4.5.1 . Wind clumping 

So far, we neglected wind clumping in our analysis. Since most 
of the nitro gen lines are formed in the intermediate or outer 
photosphere]^ at least for not too extreme winds as considered 
here, they should remain rather unaff'ected by direct clumping 
effects, though indirect effects could be important (see Paper 1). 
Clumped winds have lower mass-loss rates compared to their 
unclumped counterparts, by a factor of 1 / ^Jfc\ if the clumping 
factor /ci were radially constant, which could influence both the 
N III and the N iv emission lines, due to their sensitivity o n M. 
Given the multitude of evidence for wind-clumping (e.g., iPulsl 
l2008 l and references therein), it is necessary to examine the im- 
pact of clumping on our abundance determinations^ 

We ad opt the parametriza tion as used by IHillier & Milleil 
([T999I) and lHilher et aP (l200l . 

f(r) = + (1 - ^) exp(-v(r)/vei) (1) 

where f(r) is the volume filling factor0/t)o its asymptotic value 
(if Vco » Vci) and v^i the velocity where the volume filling factor 
reaches values close to e if /«, « 1. To allow for maximum 
effects, we set Vd - 30 kms"', close to the sonic speed for O- 
stars, and concentrated on models with /oo =0.1, corresponding 
to M reductions by a factor of roughly 0.3 (consi stent with recent 
inves tigations allowing for macro-clumping, see lSundqvist et all 
I20T1I) . 

Since most of our sample stars display thin winds, no major 
reaction due to clumping is to be expected. Indeed, a value of 
foo =0.1 did not induce any noticeable change in the spectrum 
for the bulk of the stars. 

For stars Nl 1-038, 032, B1237, Nl 1-060, 065, 066, and 068, 
mostly the N iii triplet is affected (for BI237, also N v), requiring 
0.05 to 0.1 dex more nitrogen to preserve the previous fit quality 
of these lines. Since, on the other hand, the other nitrogen lines 
remain unmodified, the inclusion of clumping did not change our 
[N] -values for these stars, but only improves or deteriorates the 
particular representation of the triplet Unes. 

Note that this might no longer be true for the N v doublet, e.g., Nll- 

031. 

-2 which is, within the standard assumption of micro-clumping and a 
void inter-clump medium, the inverse of the clumping factor /d. 



For B1253 and Sk-70° 69, clumping of the considered 
amount has a larger effect, particularly on the Nm triplet, 
NiV/14058 and the Nv doublet. Here, a clumping factor of fco 
= 0. 1 requires that [N] needs to be increased by 0. 15. . . 0.20 dex. 

Finally, for Nil -026 and Nil -03 1 , the inclusion of clumping 
affects the stellar parameters as well. Due to the lower M, a hot- 
ter temperature (by 1 to 2 kK) is needed to preserve the He i fit. 
In case of Nl 1-026 then, [N] needs to be considerably increased, 
by 0.25 to 0.3 dex. This is the only case where we encountered 
a significant effect. For Nl 1-031, on the other hand, at least the 
'cool solution' (see Sect. |5]l remained at the previous nitrogen 
abundance. 

Besides these general effects, for a few stars, Sk-66° 100, 
Sk-70° 69, Nl 1-051, Sk-66° 18 and Nl 1-065, the inclusion of 
clumping (of a lesser degree than foo = 0. 1) favours a better fit to 
Niv^4058, see Sect|5] 



4.5.2. Background abundances 

One of the central results from Paper I was that the N m triplet 
emission increases with decreasing background metallicity, Z, 
due to re duced line-blocking. In this investigation we assume, 
following iMokiem et al] (l2007bl) . a global Z of the LMC, Z = 
0.5 Zq. Since this value is somewhat c ontroversial, a nd a Z = 
0.4 Zq might be appropriate as well (e.g.. lDufouijl984 '). we need 
to test the impact of this difference on the derived nitrogen abun- 
dances. 

Overall, such lower background Z does not produce any ex- 
treme changes. As expected, the triplet emission increases, re- 
quiring a roughly 0.05 dex lower abundance to recover the pre- 
vious fits. Interestingly, we also found that NiV/14058 and the 
Nm quartet lines tend to more emission and weaker absorp- 
tion, respectively, but to a lesser extent. The effect on the re- 
maining nitrogen lines is marginal. Note that a lower [N] value 
due to lower background abundances would partly cancel with 
the corresponding increase in the derived [N] due to moderately 
clumped winds. 



5. Comments on the individual objects 

In the following, we give specific comments on the individual 
objects, regarding peculiarities and problems found during our 
analysis. We separate between B-/late O-stars and (hotter) O- 
stars, and sort by luminosity class and spectral type, starting at 
the hotter side. All nitrogen line fits (including corresponding 
limits according to Table |6]l are displayed in Appendix. |C] ex- 
cept for the objects Nl 1-072, Nl 1-032, and BI237 which have 
been included in the main paper, since they are exemplary for 
objects with different features. All spectra are coiTected for ra- 
dial velocity shifts. 

We selected those lines that were clearly visible, at least 
in most of the cases. In particular, for Teff < 35 kK, we used 
NuAA 3995, 4447, 4601, 4607, 4621, 4630 and NiiiAA 4003, 
4097, 4195, 4379, AA4634 - 4640 - 4642, and ^^4510 - 4514 - 
45 1 8. For Teff > 35 kK, we analyzed the lines from N iii/N iv/N v: 
Nm as before, NivAA 4058, 6380, and N vAAA6Q3 - 4619 (see 
Table |2]i. Additionally, the Niv multiplets around 3480 A and 
7103-7129 A have been used for the field stars observed with 
UVES (except for Sk-70° 69 where only N iv^3480 is available). 
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5.1. Late O- and B-supergiants/giants 

5.1.1. Supergiants 

Nll-029 - 09.7 lb. This is the only O-supergiant within our 
sample. Since it is of late nature, we discuss it here together with 
the B-supergiants. No major problem has been encountered, and 
the largest difference with respect to Mok07 concerns the large 
/3 - 1.63 derived in their analysis. We opted for a lower value, 
y6 = 1.23, still at the limit of prototypical values. To compensate 
for this modification, M needs to be somewhat increased. 

Figure IC.ll shows the best fit for the nitrogen lines. An 
abundance of [N] = 7.43+0.15 has been inferred, mainly from 
N 11/13995 and the Niii quartet lines. The bulk of the Nii lines 
are not helpful since there are blended by O ii (Table |2]i. This is 
also the case for most of the N iii lines. 

The most interesting feature, however, is the discrepancy for 
N III/14634H which shows an almost completely refilled pro- 
file but is predicted to be in absorption. To synthesize a pro- 
file with E.W. X would require a higher Teff or a lower logg, 
inconsistent with the He ionization equilibrium. This seems to 
be the first observational evidence for one of the problems dis- 
cussed in Paper I. For a certain temperature range (around 30 
to 35 kK), FASTWiND spectra predict too few emission in the N iii 
triplet, compared to results from cmfgen, due to (still) missing 
overlap effects between the N iii and O m resonance lines around 
374 A (which are treated consistently in cmfgen; for details, see 
Paper I). 

For further tests, we calculated a cmfgen model at the same 
stellar/wind parameters and abundances as derived in the present 
analysis. As expected, the coiTesponding synthetic profiles are 
closer to the observations (though still not as refilled as ob- 
served). To check whether the oxygen abundance plays a signif- 
icant role, two different abundances were considered, [O] - 8.66 
and 8.30 dex (solar and factor two lower). In agreement with our 
theoretical argumentation from Paper I, such a difference did not 
affect the predicted N iii triplet emission strength. 

Nll-036 - B0.5 lb. The largest difference to Mok07 is that we 
derive a lower value for log^ (by 0.2 dex) as well as a lower T^ff 
(by 500 K), by exploiting He i/He n in parallel with the N n/N iii 
ionization equilibrium. 

The nitrogen lines are well fitted, both for N n and N iii 
(Fig. lC.2l i. The only discrepancy found relates to an underpredic- 
tion of absorption strength in N n/14607, whereas the discrepancy 
at N 11/14601 is due to an On blend. Quite a large enrichment is 
found, [N] = 7.85+0.17. 

Nll-008 - B0.7 la. Again, we adapted the rather large velocity 
field exponent derived by Mok07, j6 = 1.87, to a more typical 
value of fi - 1.30. We derive v sin / = 46 kms"', which is ap- 
proximately half the value as obtained by Mok07, and compen- 
sate by invoking a v,nac = 60 km s"' . 

The best fitting abundance, [N] = 8.08+0.11 (Fig.|C3ll, has 
been obtained by the 'curve of growth method', and yields rea- 
sonable fits except for N 11/1/14447, 4621, and Nra/14003 which 
are sUghtly overpredicted. 

5.1.2. Giants 

Nll-042 - BO III. We derive a lower projected velocity than 
Mok07, V sin / = 21 kms"', as well as a lower Teff, together 

similar to N iii/i4640 which is strongly blended by O ii. 



with a corresponding decrease of logg, for consistency with the 
N ii/N III ionization equilibrium. For this cooler solution, the he- 
lium abundance needs to be lowered as well, }He= 0.08, close to 
the LMC He baseline abundance (see Sect.|6]l. 

Due to the low rotational speed, we are able to clearly in- 
spect all Nil and Niii lines (Fig. IC.4b which are fitted almost 
perfectly. The only discrepancy occurs at Niii/14097, caused by 
a coincident On line. We are able to see two strong On absorp- 
tion lines at both sides of N iii/14640, and to the left of N in/14379 
(where the former cannot be used for the diagnostics of similar 
stai-s with rapid rotation, N 1 1 -008, N 1 1 -029, N 1 1 -036, and N 1 1 - 
045). 

A low nitrogen content has been derived, [N] = 7.00+0.15, 
consistent with the low helium abundance. 

Nll-033 - BO Illn. We find a slightly cooler Teff compared to 
Mok07, and our helium line fits suggest Yue- 0.10. 

Due to its large rotation, all nitrogen lines are severely di- 
luted and almost 'vanish' from the spectrum (Fig. lC.5l l. implying 
an upper limit of [N] - 7.28 from the N n lines. 

Nll-072 - B0.2 III. This object shows a very sharp-lined spec- 
trum, with the lowest v sin / value within our sample. A consis- 
tent solution for nitrogen and helium suggests a slightly cooler 
^eif (by 1 kK) and lower log g (by 0. 1 dex). Our best fit indicates 
FHe =0.10. 

As for the similar object Nll-042, we are able to ob- 
tain a good fit from the 'curve of growth' method, with 
[N] = 7.68+0.15, and a rather low Vmic - 2.7 kms"'. Thus, and 
in contrast to Nl 1-042, this object is clearly enriched. Note that 
Nll-072 and Nll-042 belong to different associations, LH-10 
and LH-9, respectively. 

5.2. 0-stars 
5.2.1. Giants 

Nll-026 - 02 Ill(r). This star is one of the four 02 stai's in 
our sample, together with Nl 1-031, BI237, and BI253, compris- 
ing the hottest objects. Unlike Mok07, we favour a cooler solu- 
tion, Teff = 49 kK (Mok07: 53 kK), and a somewhat lower M. 
This comparably large difference for T^ff relies both on the bet- 
ter reproduction of He 1/1447 1 and on the fit to the nitrogen lines, 
with three ionization stages visible (Fig. lC.6b . Mok07 considered 
such cooler solution as well (almost included in their error bars), 
which would have improved their fit to He 1/1447 1 , but argued 
in favour of the hotter one, accounting for the global fit quality. 
By considering nitrogen now, we find support for a lower T^ff, 
since for Teff > 50 kK the N iii lines vanish from the spectrum 
(cf. BI237 and BI253). 

The derived nitrogen abundance results from a compromise 
solution, [N] - 7.80, where this value is also the lower limit. 
Moreover, we estimate quite a large uncertainty (upper limit) of 
0.4 dex, arising from our difficulties to fit all three ionization 
stages in parallel. We favour a solution that provides a good fit 
for the N iii quartet and the N iv lines, whereas a larger abun- 
dance, [N] = 8.20, is needed to fit the N v doublet. From Fig. IC.6I 
it is clear that such a large abundance (red) is in disagreement 
with the remaining nitrogen lines. Of course, we also tried a fit 
at hotter temperatures. At T^ff - 52 kK (close to the result by 
Mok07), logg = 4.1 (which is still consistent with the Balmer 
line wings) and M = 1.5-10"''MQyr"', it is possible to fit both 
Niv and N v, for quite a similar abundance, [N] = 7.75. At this 
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temperature, however, all N m lines have vanished though. To re- 
cover them we would need to increase the abundance again, also 
by 0.4 dex. Because of the poorer prediction of He 1/14471 we opt 
for the cooler solution. Since in both cases the implied nitrogen 
abundances are similar, this does not lead to severe problems for 
our further analysis, but requires larger error bars than typical. 

Nll-031 - ON2 IIW). This star raised the most severe dif- 
ficulties in our sample when trying to fit the nitrogen lines 
(Fig. IC.7b . Already its 'ON' designation indicates strong nitro- 
gen featur es in its spectra , in this case N iv/l/l4058, 3480 (for the 
latter, see IWalborn et al. 1 12004). and the Nv doublet lines. We 
were not able to consistently fit these strong features together 
with the remaining nitrogen lines. In contrast, N iv/l6380 has al- 
most the same strength as in Nl 1-026. 

If we try to reproduce the (rather weak, but clearly visi- 
ble) Hei/i4471 in pai-allel with Niii and NiV/16380, we obtain 
T^g = 47.8 kK, a bit cooler than Nl 1-026, whilst Mok07 derived 
T^eff = 45 kK, excluding T^ff values larger than 47 kK based on 
the helium ionization analysis. 

If, on the other hand, we try to fit the problematic lines, 
we need a much higher Tgff . A consistent fit for all N iv lines 
(including NiV/16380) together with those from Nv requires 
Jeff = 56 kK, logg = 4.00, and M = 2.2-lO"''M0yr-', to- 
gether with a very high abundance, [N] - 8.30. Clearly, this 
set of stellar parameters neither reproduces the N iii lines nor 
the weak He 1/14471. Nll-031 has been previously analyzed by 
[Walborn et al. (20Q1) and Iboran & Crowther (201 1) using the 
Niv/Nv lines (without discussion of Hei and Niii). The for- 
mer authors obtained quite similar parameters, T^ff - 55 kK, 



log^ = 4.00, and a somewhat lower M = 1.0 ■ lO^^'Mgyr"', 
presumably because of a clumped wind (though clumping 
has not been mentioned). At these parameters, they derived 
[N1 - 8.00+0.18, which f or our models would still be too low. 
iDoran & Crowthe^ (1201 V) only provided T^ff within their analy- 
sis, deriving Teff - 54.7 kK for this object. 

Because of the similarities with Nl 1-026 (except for the two 
strong features), the fact that Niv/I638(|3 behaves 'normally' 
and that He 1/14471 is clearly visible, the cooler solution with 
[N] = 7.83+0.15 cannot be discarded so far We checked our 
synthetic spectra by independent cmfgen calculations. The re- 
sults are quite similar, in particularly the predicted N m emission 
lines are even a bit weaker than those produced by fastwind, 
again pointing to a cooler solution. 

We also tried to attribute the problematic feature to the pres- 
ence of X-rays, by means of cmfgen calculations including typi- 
cal X-ray strengths and distribution, but almost no effect on these 
lines has been found (basically because the line forming region 
is still inside or close to the photosphere). 

Of course, it would be helpful to consider NiV/13480 as 
well, which unfortunately is not included in our dataset. A 
by-eye compa r ison w ith the corresponding profile displayed in 
IWalborn et al.l (|2004 their Fig. 1) showed that both our cooler 
and hotter solutions are not incompatible with this spectrumPI 



which turned out to be quite 'rehable' in the remaining objects, 
even for the hottest ones. 

Interestingly, our cooler solution predicts somewhat stronger ab- 
sorption than the hotter one, even though the corresponding nitrogen 
abundance is significantly lower. 
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Thus, the nature of the strong N iv/l4058 and N v doublet fea- 
tures remains open. Of course, binarity could be a plausible solu- 
tion, where a cooler component could be responsible for He i and 
N III, and a hotter one for the intense N iv and N v lines. We note 
that this is the brightest (Mv=-5.78) object of the 0-star sample 
and that other, previously thought single early O-stars, display- 
ing both strong N v lines as well as the presence of He i in their 
optical spectra such as CygOB 2-22 and HP 9 3 129 A, we re sub- 
seque ntly resolved as binaries dWalborn et alj [2002, Nel an et alj 
l2004l) . Thus, the binarity scenario for Nl 1-031 needs to be clar- 
ified in future investigations. 

Note, however, that the other ON-stars discussed by 
IWalborn et aU l2004 LH64-16 and NGC 346-3, seem to dis- 
play very similar features, though the presence of He i is not as 
clearly visible as in our spectrum for Nl 1-031. Similar, but less 
dramatic, problems are also found for those other sample stars 
where we were able to analyze N iii/N iv/N v in parallel, namely 
Nl 1-026 and Nl 1-060 (see below). While the presence of such 
discrepancies in all these objects may point towards a less likely 
binary scenario for Nl 1-031, we note that the differences in T^ff 
for the two alternative solutions reach 8,000 K in Nl 1-031, and 
remain at moderate 3,000 K for the other two objects analyzed 
here. 

We will reconsider Nl 1-031 and other ON-stars in future in- 
vestigations, to clarify the question how it is possible to have 
weak He i and N in in parallel with strong N iv and N v. For the 
remainder of this paper, however, we discuss the results forNl 1- 
03 1 in terms of both the cool and the hot solution, without pre- 
ferring either of them. 



Nll-038 - 05 Il(f^). The parameter set derived for this stars 
is quite similar to Mok07, with only slightly lower Fhb = 0.08 
and V sin / = 100 km s"' . This star displays a peculiar He 1/1447 1 
profile of triangular shape that could not be reproduced, even if 
invoking macro-turbulent broadening (as speculated by Mok07). 

We obtain [N] = 7.85+0.15, similar to the case of Nll-032, 
by means of quite a good fit to all lines from three ionization 
stages. Besides the peculiar He 1/14471 profile, Niv/16380 shows 
contamination by the DlBs at AA 6376.08, 6379.32, and the 
N v/14603 absorption is much stronger than predicted, contrasted 
to the other component. 



Since the N iv lines are weak and rather noisy, we can only state 
that our simulations are consistent with the observations. 



Nll-045 - 09 III. Our analysis is in good agreement with 
Mok07. We confirm that a low He abundance (lower than the 
estimated LMC He baseline abundance) matches the observa- 
tions. 

The only clear N-abundance indicators are the Niii quar- 
tet lines, since most other lines are weak and the spectrum is 
noisy (Fig. IC. 10b . Niii/14097 is also weakly visible, and consis- 
tent with the quartet lines. Around this T^ff, the N iii triplet turns 
from absorption to emission, and thus this object does not be- 
long to the 'f ' category. The absence of lines, the fact that the 
star does not display a fast (projected) rotation, and the very low 
helium content is consistent with the very low nitrogen abun- 
dance, [N] = 6.98+0.20. These findings suggest that this object 
is of unevolved nature. 

5.2.2. Dwarfs 

BI253 - 02 VCCf*)). This w as one of the stars that were 
used by Wal born et al.l (l2002h to define the 02 spectral type. 
A slightly hotter solution (by 1 kK, Tes - 54.8 kK) than in 
Mok07 was obtained, by using the N rv /N v ionizatio n equi- 
libr ium. This star was also analyzed by iMassev et al.l (l2005h 
and lDoran & Crowthed (1201 1). The former authors provide only 
a lower limit on Jeff and a consistently lower logg (T^ff > 
48 kK and log^ - 3.9), both in agreement with the error bars 
by Mok07, whilst IPoran & Crowthed (1201 1 ^) derive a somewhat 
cooler solution (by 2 kK) compared to our findings. 

Figure ETT] shows that Niv/16380, NiV/13480, the Niv mul- 
tiplet around 7120 A( where we reproduce the observed emis- 
sion), and the N v lines are nicely fitted, with [N] - 7.90±0.15. 
At this Teff and M, no N iii is visible in the spectrum (see also our 
discussion on Nl 1-031). The feature located around Niii/14634 
corresponds to OiV/14632. On the other hand, we were not able 
to reproduce the rather broad emission of NiV/14058 (a larger 
V sin / is inconsistent with the remaining lines). For this star, 
we compared again with a cmfgen model at similar parameters, 
in particular with the same [N]. Contrasted to our solution, the 
width of the N iv emission line could be fitted, whilst the other 
lines indicated either a hotter temperature or a higher abundance. 



Sk-66° 100 - 06 11(f) is one of the field stars within our sam- 
ple. The inspection of the H/He fits suggests no major revision 
of the values provided by Mok07. 

Figure lC9l shows the best fit to the nitrogen lines. An abun- 
dance of [N] - 8.48+0.15 was needed to obtain a consistent fit. 
This large value agrees well with the large He content found for 
this object, Fne - 0.19. However, we also found some problems 
regarding the N iv fits, except for the multiplet around 3480 A 
where the fit is perfect. Interestingly, we were able to 'cure' these 
problems by invoking a clumped wind with /oo = 0.2, with no 
significant changes in the remaining nitrogen lines. 

Nll-032 - 07 11(f). No major problems were found for this 
star. Our solution is slightly hotter than in Mok07, and we pre- 
ferred a typical /3 value for O-stars, /3 - 0.80. 

An excellent fit is obtained for this star of the (f) category 
(Fig. ISll, resulting in [N] - 7.87+0.15 for both the Niii triplet 
and the quartet lines. The fit quality of N iii/14097 is remarkable 
as well. Note that for this star N iii/l4003 has not been observed. 



BI237 - 02 V((r)). As for BI253, the nitrogen ioniza- 
tion equilibrium favours the hotter solution proposed by 
M ok07 (r^fF = 53.2 kK ) rather than the cooler limit derived 
by iMassey et alj ( 2005h. Our resu lt is also consistent with the 
work bv fPoran & CrowtheJ (1201 lb . This object is quite similar 
to B1253 with a somewhat thinner wind. 

Figure|2]displays a good fit for the nitrogen lines. Niv/14058 
does not show a broad profile as in BI253, and we are able to per- 
form an excellent fit to this line. We derive a lesser enrichment 
than forB1253, with [N] = 7.38+0.15. 

Nll-060 - 03 V((r)). With = 48 kK, which is about 2 kK 
hotter than in Mok07, we found a consistent description of H/He 
and nitrogen. As argued by Mok07, He 1/14471 is similar to Nl 1- 
031, and they derived a similar T^ff around 45 kK for both stars. 
Our findings indicate a larger value for both stars, indicating the 
internal consistency. 

Again, we encountered the problem already met for Nl 1-026 
and Nl 1-031, i.e., it is rather difficult to fit the lines from three 
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different nitrogen ionization stages in parallel, see Fig. lC.12l We 
opted for a compromise solution with [N] = 8.20+0.30, which is 
quite large but in line with the helium enrichment. Note that we 
also estimate quite a large upper limit, motivated by the follow- 
ing reasoning. 

As for Nil -026, we considered the impact of a hotter solu- 
tion. With Teff = 51 kK and logg = 4.1, we are able to fit both 
Niv and N v at a similar abundance, [N] - 8.15. To recover the 
N III triplet lines at such hot temperatures, however, requires an 
increase of [N] by roughly 0.3 dex. Thus, we are able to derive 
a quite similar nitrogen abundance for different Teff , either using 
the N m/N iv or N iv/N v ionization equilibrium. The upper limit 
results from the condition to match either N v or N iii, respec- 
tively. 

Sk-70° 69 - 05 V((f)). For this field star, we had some prob- 
lems to reconcile H/He and N at the parameters provided by 
Mok07. Our final solution is cooler by 1 kK, and we derived 
^He = 0.14 which is lower than the Yue = 0.17 estimated by 
Mok07. 

Figure IC.13I presents the best fit for all nitrogen lines, for 
[N] - 8.05+0.15. There are only two disagreements: the 'right' 
wing of N III/14634 is predicted a bit too narrow, and N iV/14058 
is predicted to be in very weak emission, contrasted to the obser- 
vations. We consider the fit as acceptable, particularly since in a 
clumped wind this weak emission becomes almost suppressed. 

Nll-051 - 05 Vn((f)). This stai- is the fastest rotator in our 
sample, and we derived a slightly higher v sin / and cooler Teff . 
As well, we used a prototypical value of (3 - 0.8 instead of the 
rather low fi — 0.6 derived by Mok07, implying also a lower M. 

Unlike Nl 1-033, where the fast rotation removes almost all 
information, this object shows the N iii triplet in emission, which 
is quite well reproduced by our model, with [N] - 7.58+0.20 
(Fig. IC.14b . At such high v sin these lines are blended with 
Cm/1/14647 - 4650 - 4652. Since carbon is not included in our 
calculations, it is not possible to predict the right wing of the 
blended profile. The fact that the carbon triplet has the same 
strength as the N m one indicat es that this star could belong to 
the n ewly defined Ofc category dWalborn et al.ll20Tol I Sofa et aP 
1201 1 ). which seems to be strongly peaked at spectral type 05 for 
all luminosity classes. 

The only problem of our fitting procedure is found for 
NiV/14058, predicted to be in slight emission and actually not 
present in the observed spectra. Contrasted to the case of Nll- 
058 and Nl 1-065, it was not possible to circumvent this discrep- 
ancy by lowering M, since the fit to Hq. becomes inacceptable 
then. By means of a clumped wind with lower M, on the other 
hand, we can fit both H„ and obtain a better result for N iV/14058, 
whilst not compromising the remaining nitrogen lines. 

Nll-058 - 05.5 V((f)). To find a consistent solution for all H, 
He, and N lines, a very low M - 0.01-10"^ Moyr"' is required 
(similar to the case of Nll-065), and also a low logg (0.14 dex 
lower than Mok07), which is still consistent with the Balmer line 
wings but appears to be rather low for a dwarf. 

This was the only O-star that could be analyzed by the 
'curve growth method' with respect to [N] and Vmic. The de- 
rived value of v'mic = 6±3 kms 'is lower than for the other O- 
stars, though such a difference does not drastically affect the 
derived abundance, as argued in Sect l4.4l Quite a large abun- 



dance, [N] = 8.09+0.15, was determined which fits all the lines 

(Fig. EH). 

Sk-66° 18 - 06 V((f))- Our estimates agree well with those 
from Mok07. As for Sk-66° 100, we found a very large ni- 
trogen abundance, [N] - 8.48+0.15 (Fig. IC.I6I 1. Note that the 
NiV/14058 line is predicted in weak emission but appears in ab- 
sorption. Again, lowering M was not sufficient to cure this prob- 
lem. A somewhat better fit to the this line was obtained for a 
weakly clumped wind with /d = 2.3, included in the figure. 

Nll-065 - 06.5 V((f)). Our best fit to the He lines indicates 
Fee = 0.13, lower than the value derived by Mok07, Fne = 0.17. 
As already discussed in Sect. 14.31 a satisfactory reproduction of 
NiV/14058 requires a very low M - O.O5-lO"^M0yr"'. Since this 
line clearly appears in absorption, whereas our model with M 
from Mok07 predicts much too less absorption, we lowered M, 
but were able to preserve the fit to H„ and the remaining nitrogen 
lines. By the inclusion of clumping we obtained an even better 
fit quality, with [N] = 8.17+0.15 from N iii and N iv. 

Nll-066 - 07 V((f)). Teff and log g as derived by Mok07 turned 
out to be inconsistent with the N iii/N iv ionization equilibrium. 
We had considerable problems to fit both N iii and N iv lines at 
the same abundance, and particularly to reproduce the absorp- 
tion within Niv/I4058and the Niii quartet lines. To this end, a 
lower Teff was mandatory, but only rough estimates using our 
coarse grid could be obtained, resulting in Teff = 37 kK and 
logg - 3.7 dex, which are 2.3 kK and 0.17 dex lower than the 
values provided by Mok07, respectively. The gravity is some- 
what low for a dwarf but still inside the error bars assigned by 
Mo k07. W ith these values, we determined [N] = 8.17+0.20, see 
Fig. lC.18l which seems rather large for Fhc = 0.1. 

Nll-068 - 07 V((f)). As for Nl 1-066, we were only able to fit 
the H/He and N lines by means of the coarse grid. The differ- 
ences to the results by Mok07 are significant, but still consistent 
with the observations and identical to those of Nll-066 which 
has the same spectral type. 

In contrast to Nl 1-066, however, particularly the N iii triplet 
shows weaker emission, and thus a lower [N] = 7.85+0.20 has 
been found, see Fig. lC.19l 

Nll-061 - 09 V. Contrasted to Mok07, we chose a prototypical 
value for the velocity field exponent, y6 = 0.8, together with a 
larger value for M. All H/He lines could be reproduced without 
difficulties. 

The nitrogen analysis is quite similar to the case of Nl 1-045, 
and we derived [N] = 7.18+0.15 fromNiii alone (Fig. lOOb . 

Nll-123 - 09.5 V((f))- Our parameters are very similar to 
Mok07, and Fig. IC.21I shows our solution for Nii/Niii. Even 
though the rotation is not extreme, almost no nitrogen is visible, 
and the features overlapping with the 'non-existent' N ii/l4630 
and N iii/l4640 are blends by O iii and Si iii, respectively. Thus, 
we infer a very low nitrogen content, [N] = 7.00+0.15, roughly 
corresponding to the LMC baseline abundance. 

Nll-087 - 09.5 Vn. Also for this rapid rotator, we found good 
agreement with Mok07. Due to rotation, all nitrogen lines are 
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diluted into the continuum (Fig. IC.22b . and only an upper limit 
for [N] could be estimated, [N] = 7.38. 

6. Discussion 

6.1. Comparison with resuits from MokOJ 

In our discussion of the derived results, let us first concentrate 
on a brief comparison with the findings by Mok07 (see Table|5]l. 
Except for few cases, we derive somewhat cooler T^ff. To a cer- 
tain extent, this might be attributed to the improved temperature 
structure in the new fastwind version, and also to the possibil- 
ity to exploit the nitrogen ionization equilibrium. The largest 
changes concern Nl 1-026 (4.3 kK cooler), and Nl 1-031 (al- 
ready the cooler solution is 2.8 kK hotter). For the earliest stars 
in our sample (two 02 dwarfs and two 02 giants), we determine 
a temperature range of 47.8 < T^g < 54.8, quite similar to 
Mok07, but now inc luding Nl 1-031, sin ce we infer hotter solu- 
tions for this object. iMassev et aP (l2005h found a similar range, 
using two dwarfs and two giants. Gravities changed in parallel 
for most of the cases, with 0. 10 and -0.43 dex as largest positive 
and negative difference, respectively. Stellar radii show a very 
good agreement, even for the two stars with the most extreme 
changes in Teff (less than 5% difference in /?»). 

Regarding M, we derive values which are typically lower by 
less than a factor of two (A log M « -0. 1 . . . -0.4 dex), with a 
maximum change of - 1 .2 dex for Nil -058 based on our analysis 
of the nitrogen lines. The agreement of the resulting Fhc values 
is good, except for two stars with differences considerably larger 
than the adopted errors. For both stars (Sk-70° 69 and Nl 1-065), 
we find a lower helium content. 

The largest differences relate to v sin / and Vmic. Differences 
around 30-40% in v sin / stress the importance of obtaining this 
parameter in a separate step, when using an automated fitting 
method. The substantial differences in Vmic, on the other hand, 
should not be regarded as worrisome though. To a major part, 
Vmic has not been literally derived during this work, but was only 
adopted (as Vj^ic- 10 kms"'), where the resulting fit quality did 
not indicate any problems with this value, within + 5 kms"'. 
Only for four mostly cooler stars we were actually able to infer 
more robust estimates, indicating quite a low v^ic ~ 5 kms"'. 
Note, however, that the latter value refers to nitrogen lines only, 
and inconsistencies in Vj^ic from H/He (used by Mok07) and 
metal lines have been found already in various studies. 

Since the topic of a potential relation between Vmic and stel- 
lar type (log^!) is of recent interest!^ because it might indicate 
(together with other evidence) the presence of sub-surface con- 
vection (iCantiello et al.ll2009h . a more thorough investigation is 
certainly required. A derivation from light elements will become 
difficult for the hotter O-stars though, due to the restricted num- 
ber of visible lines and the complex formation mechanism of the 
ubiquitous (photospheric) emission lines. Here, it will become 
advantageous to exploit the i nformation contain ed in the numer- 
ous UV Fe and Ni Hnes (e.g. jHaser etalJll998h . 

6.2. Overiap witli B-star nitrogen anaiyses 

To ensure the consistency between O-star nitrogen abundances 
from this (and upcoming) work and previous results from B-stars 



e.g., M okO? found a weak correlation for objects with log^ < 
3.6; se e also lKilian et"an dl99lh . rGies & Lambert (ll992h.lDaflon et alJ 
( 120041) for similar results for Galactic B-stars and lHunter et all ilOOlh 
for LMC B-stars. 



Table 7. Comparison bet ween nitrogen abun dances derived by 
means of fastwind and bv lHunter et all (l2009l) . using their stellar 
parameters and Nii diagnostics alone. 



Star 


Sp. Type 


[N] t.w. 


[N] Hunter 


Nil -029 


09.7 lb 


7.20 


7.10 


Nl 1-036 


B0.5 lb 


7.86 


7.76 


N 11-008 


B0.7 lb 


7.89 


7.84 


Nl 1-042 


BO III 


6.99 


6.92 


Nl 1-072 


B0.2 III 


7.43 


7.38 



(using different codes, model atoms and analysis methods), a 
more thorough inspection of the cooler objects is certainly re- 
quired. Indeed, we are able to compare with alternative nitro- 
gen abundarices fr om some overlapping objects (compiled by 
Hun ter et al.l2009t see Table|6]l, which base on stellar parameters 
obtained by means of tlusty an d using the Si m/Si iv ionization 
equiUbrium dHunter et alJl2007l) . 

Already in the latter work. [Hunter et al] (l2007h realized that 
their Teff-values were somewhat lower than corresponding re- 
sults from Mok07 (who used exclusively H and He), but also 
pointed out that Teff estimates based on Heii/14541 would be in 
much closer agreement. This is even more true regarding our 
'new' values, which lie in between the Hunter et al. estimates 
from Si and those from Mok07. Such differences in T^g derived 
either from metals or from H/He are somewhat disturbing, due 
to their influence on the metallic abundances when using lines 
from one ion ization stage only. 

Anyhow, iHunter et al.l (l2007l) decided to keep their cooler 
solution, to preserve the internal consistency of their analysis. 
Consequently, the nitrogen abundances derived in the present 
study are systematically larger, due to our higher T^ff (leading 
to intrinsically weaker N ii lines), with AFetf ~ IkK for Nl 1-008 
and Nl 1-072, and ^ 2kK for Nl 1-036. For Nl 1-042, the T^ff 
is rather similar, and for this object the derived [N] -values in- 
deed overlap by better than 0. 1 dex. Also for Nl 1-029, our Teff 
is rather similar, but we derive a 0. 1 dex lower log g, which leads 
to weaker N ii lines. 

To check our model atom and our diagnostic approach, we 
performed an additional analy sis by reproducing the conditions 
adopted bv lHunter et al.l(l2007l) . i.e., we used their stellar param- 
eters (with negligible M to mimic tlusty models) and N ii lines 
only, within the 'curve of growth' method and adopting their N ii 
equivalent widths and uncertainties. For Nl 1-029 and Nl 1-042, 
we could only perform a 'by eye' estimate because Hunter et al. 
considered one Nii line alone. Corresponding results are com- 
pared in Table |7] and the agreement is almost excellent. 

Thus, we conclude that the consistency of nitrogen abun- 
dances in the overlapping B- and O-star regime is satisfactory, 
and that the different codes and methods produce a reasonable 
agreement. However, there is a slight offset on the order of 0. 1 
to 0.2 dex, which we attribute mostly to different effective tem- 
peratures. Since our analysis is based on both N ii and N iii lines 
(in contrast to Hunter et al.), and has been performed in parallel 
with the analysis of H/He, we prefer our values though. 

6.3. Nitrogen abundances 

Figure |8] summarizes the basic outcome of our analysis, by dis- 
playing the derived nitrogen abundances as a function of spectral 
type and he lium content, t ogeth er with the LMC baseline abun- 
dance from iHunter et al.l ( l2007h . Evidently, there are only few 
cooler objects located close to the baseline, whereas the major- 
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Fig. 8. Nitrogen abundances derived for our LMC sample. Upper 
panel: As a function of O-star spectral type ('10-12' correspond 
to B0-B2; for Tsg, see Table |6]l. Lower panel: As a function of 
helium content, 12 + log(He/H) = 12 + log(FHe)- Red: luminos- 
ity class I-II; blue: III-V. Arrows indicate upper or lower limits. 
Alternative solutions for Nl 1-031 (see Sect.|5]l are indicated by 
a triangle and an asterisk for the cooler and hotter solution, re- 
spectively. The estimated (1-cr) error for lo g(YHe) is 0.05 dex . 
The dotted lines indicate the LMC nitrogen ( Hunter et al.ll2007h 
and the average helium (see text) baseline abundances. Some 
objects have been slightly shifted horizontally, because of better 
visibility. 



ity of the objects (independent of luminosity class!) is strongly 
enriched, with [N] in between 7.5 and S.lFlFive objects display 
extreme enrichment, with [N] from 8.17 to 8.5, which is close 
to the maximum nitrogen content given by the CNO equilibrium 
value for nitrogen, [N]n,ax ~ 8.5. However, this is well above the 
enrichment reached fo r a 40 Mp star w ith an initial rotational 
velocity of 275 km s ' dBrott et"d]|201 1a'). 

The lower panel of Fig. [8] is more promising though. 
There seems to be a strong correlation between the ni- 
trogen and the helium enrichment, here displayed loga- 
rithmically. The LMC helium abundance should be lo- 
cated, in terms of number fraction, around Yue = 0.08- 
0.094 con-esponding to [He] = 10.90-10.97 (Rus sell & Dopital 



Note that most stars have 'normal' or only moderately enriched he- 
lium abundances (see below) so that the high [N] cannot be an effect of 
decreasing H content. 
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Fig. 9. 'Hunter-plot' displaying the nitrogen abundan ce vs. pro- 
j ected r otational speed. Population synthesis from iBrott et al.l 
( 1201 Ibh . for Teff > 29kK and a magnitude limit V < 15.29, 
shown as a density plot in the background. The color coding cor- 
responds to the number of stars per bin, with binsize 5 kms 'x 
0.04 dex. Overplotted are data from this study. Black: Nil stars; 
blue: field stars. Circles, diamonds and squares correspond to 
objects with low, intermediate and strong helium enrichment, 
respectively (see text). Alternative solutions for Nl 1-031 (inter- 
mediate He enrichment) as in Fig.[8l 



19901: iMaeder & Me vnetl 120011: IVermeii & van der HulstI 120021 



PeimbertH200'3r iTsamis et al.ll2003h . and agrees quite well with 
our minimum values for the derived helium abundance. We 
found only five stars with considerable enriched nitrogen close 
to this value, three (super-)giants and two dwarfs, but note also 
the attributed uncertainty in helium content. Except for these ob- 
jects, the correlation is almost perfect, and there is a certain clus- 
tering around the pair [He] = l 1.0/[N]=8.0. 

A somewhat different view is provided in Fig. |9] which dis- 
plays the so-called 'Hunter-plot', nitrogen-abundance vs. pro- 
jected rotational speed, for all our sample stars with Teff > 29 kK. 
N-1 1 stars are indicated in black, field stars in blue. Circles, di- 
amonds and squares correspond to objects with low (Fee < 0.1), 
intermediate (Yue - 0.1), and strong (Yue > 0.1) helium enrich- 
ment, respectively. 

The background of this fig ure consists of resu lts from the 
recent population synthesis by iBrott et aH (1201 Ibl) . for all ob- 
jects with Teff > 29 kK, and a magnitude limit (corresponding 
to our sample) of V < 15.29, shown as a density plot. The un- 
derlying simulation assumes a rather broad Gaussian rotational 
velocity distribution as derived for LMC early-type stars, peak- 
ing at 100 kins" 'wi th a standard deviation of cr = 140 kms ' 
(iHunter et al.ll2 008b. 2009)0 and random inclinations. 

Such diagrams ([N] vs. v sin /, compared with evolutionary 
calculat ions) have been presented the first time bv lHunter et al.l 
(l2008ah . to summarize the outcome of the B-star analyses within 
the FLAMES-I survey, and to investigate the predicted effects 
of rotational mixing. One of their major findings was the un- 
expected presence of a significant number of objects with slow 
rotation and large enrichment, not predicted by (single-star) the- 
ory, so-called 'group 2' objects. 

In the O-star case now, this problem becomes even more se- 
vere. We refrain here from a detailed statistical analysis, since 



Supergiants and stars above 25 Mq have been discarded from their 
analysis, to avoid effects from mass loss induced spin-down. 
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the number of investigated objects is too low, and postpone this 
objective until the results from the FLAMES Tarantula survey 
(with more than 200 'useful' O-stars) have become available. 

Nevertheless, the trend is obvious. Roughly one third of the 
objects are located at positions where they should be expected 
(those at the baseline and the 'diagonal'), another one-third is lo- 
cated at the predicted upper limit, and the last one-third (beyond 
[N] - 8.0) extends to very large values where the predicted popu- 
lation density is almost zero. Let us note that the two objects with 
the largest [N] enrichment 8.5 dex, which is (incidentally?) 
just the maximum nitrogen content given by the CNO equilib- 
rium value) are two field stars, Sk-66° 100 and Sk-66° 18, one 
06 giant and another 06 dwarf. Both stars did not present any 
difficulties in the nitrogen analysis, thus indicating a reasonable 
quality. 

In terms of the original Hunter et al. sample, roughly two- 
third of our objects would be denoted by 'group 2' Pi The cor- 
responding number of objects is so large that inclination ef- 
fects w.r.t. V sin ; should be irrelevant. Interestingly, however, 
the corresponding He-abundances are in line with our findings. 
The first group has a low abundance (Fig. |9] circles), the sec- 
ond group mostly consists of enriched objects (diamonds), and 
the third one comprises objects with considerable He enrichment 
(squares). Thus, in parallel with the derived correlation between 
observed nitrogen and helium content, the discrepancy between 
observations and theory becomes the stronger the larger the He- 
abundance is. 

In evolutionary models is the amount of He transported to 
the surface strongly controlled by the parameter ^, which de- 
scribes the inhibiting effect of mean molecular weight gradients 
(in t his case, the H-He gradient) on the tr ansport of elements 
(see lHeger&Lang er 2000). The models of lBrott et all (1201 lal) 
have a dopted ^ =0.1, from an earlier calibration of Yoon et aTl 
(l2006l) . Lowering the sensitivity to the mean molecular weight 
barrier would increase mixing of both, nitrogen and helium, to 
the surface, but also reduce the minimum mass and velocity re- 
quired for chemical h omogeneous evolution in the models (see 
also the discussion in^Heg er & Lange3l2000 l). Given the present 
values of [N] and Yue, it might be possible to derive further con- 
straints on in future work. 

7. Summary 

In this paper, we investigated the NiV/14058 emission line for- 
mation, determined the nitrogen abundance of a sample of 25 
LMC O- and early B-stars, and performed a first comparison 
with corresponding predictions from stellar evolution including 
rotational mixing. The results of this work can be summarized 
as follows. 

1. In O-stars, the dominating process responsible for the Niv 
line emission is the strong depopulation of the lower level 
by the 'two-electron' transitions 3p — > 2p^, of (mostly) pho- 
tospheric origin. This drain increases as a function of M, 
because of increasing ionizing fluxes (which are coupled to 
the Hen continuum), leading to more depopulation of the 
ground and the coupled 2p^ states. Resonance lines (as for 
the N m emission triplet) do not play a role for typical O-star 
mass-loss rates and below. 

Since in addition to nitrogen there are many other elements 
which display optical line emission in the hot star regime 



In contrast to the B-star 'group 2' objects, however, the devia- 
tions between predictions and 'observations' for some of the objects 
are much more extreme. 



(C, O, Si), it might be suspected that similar processes 
might be invoked, because of similar electronic configura- 
tions/transitions. 

2. To infer the nitrogen abundances, we re-determined the stel- 
lar and wind parameters, by means of 'by eye' fits, starting 
with the values provided by Mok07, but exploiting in par- 
allel the nitrogen ionization equilibrium and deriving v sin / 
in a first, separate step. Moreover, we accounted for extra 
line-broadening expressed in terms of Vmac- In addition to 
systematically lower v sin /, we also derived mostly lower 
Teff (partly because of using an improved fastwind version) 
and thus log g, but differences to Mok07 are generally small, 
except for few objects. 

3. Based on these parameters, we derived nitrogen abundances, 
mostly by varying the abundance and comparing with all ni- 
trogen lines present in the spectrum. In a few cases, we were 
able to estimate [N] and v,nic in parallel, by means of a curve- 
of-growth method. 

4. Again in most cases, we found no problems in fitting the 
nitrogen lines, and reproduced the 'f features quite well. 
Only for some of the (hotter) objects where lines from all 
three stages, N m, N iv and N v, are visible, we needed to aim 
at a compromise solution. Real problems were encountered 
for one star, Nl 1-031 (ON2 Ill(r)), where only either Hei, 
N m and N iv/l6380 (at cooler Teff) or N iv and N v (at higher 
Tetf) could be fitted in parallel. The difference in the derived 
Tetf amounts to 8,000 K, which is far from satisfactory, and 
requires future effort to resolve the problem. A solution in 
terms of binarity, though somewhat unlikely, cannot be ruled 
out so far. 

5. For some cooler objects already analyzed bv iHunter et al.l 
poo?) by means of tlusty using Nii lines alone, we found 
differences in [N] on the order of 0.1 to 0.3 dex, with larger 
values from our analysis. These differences could be exclu- 
sively attributed to diff'erent stellar parameters, mostly T^ff. 
Overall, however, are the nitrogen abundances in the over- 
lapping B- and O-star domain consistent within a reasonable 
error. 

6. Within our sample, we found only three cooler objects close 
to the LMC nitrogen baseline abundance, [N]baseiine - 6.9. 
The majority of the analyzed O-stars (independent of lu- 
minosity class) seems to be strongly enriched, with [N] = 
7.5 to 8.1. Five objects indicate an extreme enrichment, with 
[N] = 8.17 to 8.5. 

7. There is a rather good correlation between the derived nitro- 
gen and helium surface abundances. 

8. Comparing the nitrogen abundances as a function of v sin / 
with tailored evolutionary calculations, we found a signifi- 
cant number of highly enriched, low v sin / ('group 2') ob- 
jects. Interestingly, the correlation between He and N be- 
comes also visible in this comparison: Whilst most objects 
with unenriched He are located just in the region where the 
predicted population density is largest (accounting for selec- 
tion effects), objects with enriched He are located at the up- 
per limit of this distribution and above, and particularly those 
with the largest He enrichment lie well above this limit. 

Due to the low initial (baseline) nitrogen abundance, the detec- 
tion of strong nitrogen enrichment in the bulk of O-stars might 
indicate that efficient mixing takes place already during the very 
early phases of stellar evolution of LMC O-stars. Nevertheless, 
it would be premature to draw firm conclusions from our re- 
sults, since the sample size is still low. Upcoming results from 
the VLT-FLAMES Tarantula survey (which will be derived in 
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a similar way as presented here, drawing from our experience) 
will enable a more complete view. In particular, the determina- 
tion of O-star nitrogen abundances in the LMC will place very 
tight constraints on the early evolutionary phases of O-stars and 
thus on the theory of massive star evolution. 
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Appendix A: Details of the nitrogen model atom 

This section provides some details of our N ii, N iv and N v 
model ions (corresponding material for Nm has been akeady 
presented in Paper I). Configurations and term designations are 
outlined in Tables O] and |A3] whilst Figs. lOllAiH and 
IA.3l displav the Grotrian diagrams for the Nii singlet and triplet 
system (the quintet system comprises five levels only), the N iv 
singlet and triplet system, and the N v doublet system, respec- 
tively. In these figures, important optical transitions as given in 
Table|2]are indicated as well. 
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Fig.A.l. Grotrian diagrams for the Nii singlet (upper panel) 
and triplet (lower panel) systems. Level designations refer to 
Table lA.ll Important optical transitions are indicated by green 
lines and numbers referring to entries in Table 12 



Table A.l. Electronic configurations and term designations of 
our N II model ion. The level numbers correspond to the entries 
in the Grotrian diagrams in Fig. lA.ll for the singlet and the triplet 
terms. 
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Table A.l. Electronic configurations and term designations of 
our N IV model ion. The level numbers correspond to the entries 
in the Grotrian diagrams in Fig. lA.2l 
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Fig. A.2. As Fig. lA.ll but for the N iv singlet (left) and triplet (right) system. Level numbers refer to Table I A. 2] 
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Fig. A.3. As Fig. lA.ll but for the N v model ion. Level numbers 
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Appendix B: Tests of the Nii model ion 

Figures IB.ll to IB. 7! refer to tests of our N ii model ion, as de- 
scribed in Sect. 12.2.11 Figure IB.ll compares electron tempera- 
tures and densities for B-star parameters calculated by fastwind 
and TLUSTY, whilst Figs. lB.2 lto lB.7l compare corresponding syn- 
thetic N II line profiles from these two codes and from calcula- 
tions by Przybilla et al. (priv. comm.). 



Table A.3. Electronic configurations and term designations of 
our N v model ion. The level numbers correspond to the entries 
in the Grotrian diagram Fig. lA.3l 
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Appendix C: Line fits for individual objects 

Figures. IC. ll to lC.2^ displav the observed (green) and best-fitting 
optical nitrogen spectra (black) for all our objects, except for 
Nil -072, Nl 1-032, and BI237 which are contained in the main 
section. Blue and red spectra show corresponding synthetic line 
profiles with [N] at the lower and upper limit, respectively. For 
Nil -031 (Fig. IC.7b . we show the fits corresponding to the two 
alternative solutions for this star (see Sect. |5]l. For details on the 
line fits, see Sect. |5] and for adopted stellar parameters and de- 
rived nitrogen abundances inspect Table |4] and Table |6] respec- 
tively. All fits are based on unclumped winds except explicitly 
stated. 
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Fig.B.l. FASTWiND (black) vs. tlusty (red): comparison of electron temperature and electron density as a function of Rosseland 
optical depth, for the six models corresponding to Figs. lB.2l to lB.7l 
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Fig.B.2. Comparison of important optical Nii line profiles for a model with Teff = 20 kK and \ogg - 3.0, for models 
- FW (black), FW2 (red), TL (green) and Prz (blue). Note that N ii ^ 3995 is not present in the BSTAR2006 grid. 
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Fig.B.3. As Fig. 121 but for T^s = 20 kK and \ogg = 2.5, for models FW (black), FW2 (red), and TL(green). 



J.G. Rivero Gonzalez et al.: Surface nitrogen abundances for O-stars in the Large Magellanic Cloud, Online Material p 5 




3994.0 3994.5 3995.0 3995.5 3996.0 
wavelength (A) 
NII4621 




4620.0 4620.5 4621 .0 4621 .5 4622.0 4622.5 
wavelength (A) 
NII5005 




5004.0 5004,5 5005.0 5005.5 5006.0 
wavelength (A) 




1.1 
1.0 
0.9 
0.8 

0.7 
0.6 

4446.0 4446.5 4447.0 4447.5 4448.0 
wavelength (A) 
NII4630 




4629.5 4630.0 4630.5 4631.0 4631.5 
wavelength (A) 
NII5679 




5678.5 5679.0 5679.5 5680.0 5680.5 
wavelength (A) 



o 


1.1 i 

1.0 r 


o 




Q. 


0.9 r 


c 








0.8 r 


E 




Qi 


0.7 r 




06 i 




4600.5 4601.0 4601.5 4602.0 4602.5 
wavelength (A) 
NII4643 




4642.0 4642.5 4643.0 4643.5 4644.0 
wavelength (A) 
NII6482 




6481.0 6481.5 6482.0 6482.5 6483.0 
wavelength (A) 



Fig.B.4. As Fig.lEll but for Teif = 24 kK and \ogg^ 3.0. 
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Fig. B.5. As Fig. 13] but for T^a = 24 kK and log g = 2.75. 
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Fig.B.6. As Fig.lJl but for T^f^ = 27.5 kK and Xogg = 3.5. 
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Fig.B.7. As Fig. 13] but for Teif = 27.5 kK and \ogg = 3.0. 
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Fig. C.16. Sk-66° 18 - 06 V((f)). Magenta spectra correspond to a weakly clumped model. For details, see Sect.|5] 
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Fig. C.18. Nll-066 - 07 V((f)) 
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Fig, C.19. Nl 1-068 - 07 V((f)) 
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